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ABSTRACT 

By using updated stellar low mass stars models, we can systematically inves- 
tigate the nucleosynthesis processes occurring in AGB stars, when these objects 
experience recurrent thermal pulses and third dredge-up episodes. In this pa- 
per we present the database dedicated to the nucleosynthesis of AGB stars: the 
FRUITY (FRANEC Repository of Updated Isotopic Tables & Yields) database. 
An interactive web-based interface allows users to freely download the full (from H 
to Bi) isotopic composition, as it changes after each third dredge-up episode and 
the stellar yields the models produce. A first set of AGB models, having masses 
in the range 1.5 < M/Mq < 3.0 and metallicities 1 x lO'^ < Z <2x 10"^ 
is discussed here. For each model, a detailed description of the physical and 
the chemical evolution is provided. In particular, we illustrate the details of the 
s-process and we evaluate the theoretical uncertainties due to the parametriza- 
tion adopted to model convection and mass loss. The resulting nucleosynthesis 
scenario is checked by comparing the theoretical [hs/ls] and [Pb/hs] ratios to 
those obtained from the available abundance analysis of 5-enhanced stars. On 
the average, the variation with the metallicity of these spectroscopic indexes is 
well reproduced by theoretical models, although the predicted spread at a given 
metallicity is substantially smaller than the observed one. Possible explana- 
tions for such a difference are briefly discussed. An independent check of the 
third dredge-up efficiency is provided by the C-stars luminosity function. Con- 
sequently, theoretical C-stars luminosity functions for the Galactic disk and the 
Magellanic Clouds have been derived. We generally find a good agreement with 
observations. 

Subject headings: stars: AGB and post-AGB — physical data and processes: nuclear 



-4 - 



reactions, nucleosynthesis, abundances 
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Introduction 



Asymptotic Giant Branch (AGB) stars are among the most efficient polluters of the 
Interstellar Medium (ISM). Due to the simultaneous action of internal nucleosynthesis 
(p, a and n captures), convective mixing and a strong stellar wind, these stars provide a 
substantial contribution to the chemical evolution of the Universe. Light (e.g. He, Li, F, C, 
N, Ne, Na, Mg) and heavy (from Sr to Pb) elements can be produced. In particular, the 
activation of neutron sources in their He-rich zone makes AGB stars a special site for the 
stellar nucleosynthesis. The main neutron source is the ^^C(a,n)^^0, which is particula rly 



impo rtant in low mass AGB stars {M < 3Mf;^) (IStraniero et al 



1998 



for a review, see 



AGB St 


ars (j 


Herwie 


2005 



and 



Straniero et al. 



1995 



Gallino et al. 



20061 ). During the Thermally 



Pulsing AGB phase (TP-AGB), a recurrent penetration of the convective envelope in the 
H-exhausted core (Third Dredge- Up, hereinafter TDU), brings a few protons in a region 
where the chemical pattern is dominated by helium (80%) and carbon (20%). Then, when 
the temperature rises in this region, a ^'^C pocket forms where the ^'^C(Q;,n)^^0 reaction 
takes place. Such a ^'^C burning gives rise to the so-called main component of the s-process 
nucleosynthesis. It occurs in radiative conditions on a relatively large timescale (some 10^ 
yr), and produces a rather low neutron density (10^-10^ neutrons/cm^). A second neutron 



source may be a ctivated within t 



22 



Ne 



a,n) Mg (iCameron 



1960 



le convective zone generated by a Thermal Pulse (TP), the 



Truran fc Iben 



1977 



Cosner et al. 



198G 



Iben fc Renzini 



19831 ). In this case, the timescale of the neutron capture nucleosynthesis is definitely shorter 
(some yr) and the average neutron density is higher (10^-10^^ neutrons/cm'^). This second 
neutron burst is only marginally activated in low mass stars of solar-like metallicity, while 
it is expected to be more efficient in more massive AGB stars and a t the low metallicities 



typically found in the Gala c tic halo or in dwarf sp heroidal galaxies (IStraniero et al. 



Karakas fc Lattanzio 



2006 



Cristallo et al. 



20091). 



2000; 
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From iLambert et al. I (Il995l ). the theoretical scenario here recalled has been largely 
confirmed by the abundance analysis of AGB stars undergoing TDUs. Theoretical stellar 
models are essential to interpret the continuously growing number of spectroscopic data 
collected with the latest gener ation of modern telescopes, from very metal-poor C-rich 



1159 stars) (jWerner fc Herwig 



2006 



Werner et al. 



of measured isotopic ratios in pre-solar SiC grains (jHoppe et al. 



Zinner et al. 



2006 



Barzyk et al. 



2007 



1 

fN-tvpe) stars (lAbia et al. 


2001 


. 2002 


. 2008. 2010: 


de Lavernv et al. 2006: Zamora et al. 


2009). Ba-stars (Allen k Barbu\ 


' 2006 


: Smilianic et 


al. 


2007: Liu et al. 


2009 


, post-AGB 
), planetary 


stars (van Winckel & Revniers 


2000; ] 


■leddv et al. 


2002; 


Revniers et al. 


2004 


nebulae ( 


SharDee et al. 


200 


7; S1 


;erline 


& Dinerstein 


200 


3) and born-again AGB stars (PC 



2009h. Simi 



Davis et al. 



arly, the growing number 



1997 



Lugaro et al. 



1999 



20091 ) require specific calculations to 



interpret them. Nucleosynthesis predict i ons are also needed t o built up chemica 



models (see e.g. 



Travaglio et al. 



1999 



Romano et al. 



2010 



Kobayashi et al. 



evolution 



2010|) 



All these studies require updated theoretical predictions of the nucleosynthesis occurring 
in AGB stars belonging to very different stellar populations, with different masses and 
initial compositions. For these reasons, we decided to create a complete and homogeneous 
database of AGB nucleosynthesis predictions and yields. 

During the last 20 years, our working group has produced several improvements in the 
physical description of AGB stellar interiors and has developed the numerical algorithms 
needed to follow the complex interplay between n uclear burning, convec tion and mass 



loss. All these efforts converged in a recent paper (ICristallo et al. 



2009 



, hereafter Paper 



I), where stellar models calculations of 2 Mq AGB stars of various metallicities were 
presented. These calculations were based on a full nuclear network, from H to Bi, including 
700 isotopes and about 1000 nuclear processes among strong and weak interactions. In 
this paper, we present the FRUITY database (FRANEC Repository of Updated Isotopic 
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Tables & Yields), which is available on the web pages of the Teramo Observatory (INAfJ^. 
This database has been organized under a relational model through the MySQL Database 
Management System. This software links input data to logical indexes, optimizing their 
arrangement and speeding up the response time to the user query. Its web interface has 
been developed through a set of Pero scripts, which allow to submit the query strings 
resulting from filling out appropriate fields to the managing system. It contains our 
predictions for the surface composition of AGB stars undergoing TDU episodes and the 
stellar yields they produce. Tables for AGB models having initial masses 1.5 < M/Mq < 
3.0 and 1 x 10^'^ < Z < 2 x 10~^ are available. FRUITY will be expanded soon by including 
AGB models with larger initial mass and/or lower Z. 

In Sections [2] and [3] of the present paper we describe the stellar models and the related 
nucleosynthesis results. In Section H] we address the main uncertainties affecting our models 
while comparisons with available photometric and spectroscopic data are discussed in 
Section [61 Conclusions are drawn in Section [71 



The FRANEC code 



The stellar models of the FRUITY database have been obt ained by means of the 



Chiefs et al. 



FRANEC code (Frascati RAphson-Newton Evolutionary Code - 
most important upgrades nee ded to follow the AGB e volution and nu c. 
been extensively illustrated in 



19981). The 



Straniero et al. 



fl2006[ ) 



Crist alio et al. 



eosyn thesis have 



(l2007f ) and in Paper 



I. 



We have computed 28 evolutionary models, for different combination of masses (1.5, 



-"^ http : / / www. oa-ter amo. inaf . it / fruity 

^Perl is a general-purpose, object-oriented programming language. 
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2.0, 2.5 and 3.0 M©) and metallicities [Z = 1.0 x 10"=^ with r=0.245; Z = 3.0 x 10'^ 
and Z = 6.0 X lO^^ with r=0.260; Z = 8.0 x lO^^ and Z = 1.0 x 10"^ with r=0.265; 
Z = l.A X 10"^ = Zq and Z = 2.0 x 10"^ with F=0.269). In the present paper only scaled 
solar compositions are considered, whose relative distribution has been derived according to 



the solar compilation by 



Lodders 



(120031). 



ated for a scaled solar 



Cristallo et al. 



(120071 ). to 



In addition to the commonly adopted opacity tables calcu_ 
composition, we use C-enhanced opacity tables, as described in 
take into account the effects of TDU episodes. The most relevant atomic and molecular 
species have been considered. This improved radiative opacity substantially affects the 
stellar radius and the mass loss rate. This in turn, shor tens the AGB l ifetime hence reduci ng 



the amount of material returned to the ISM (see also 



Marigo 



2002 



2009 



Ventura fc Marigo 



Weiss fc Ferguson 



20091 ) . Our code inte rpolates between 



C-en hancement s [ft) depend on the metallicity ( 1 Crist alio et al 



opacity tables whose 



2007 



Lederer &: Aringer 



20091 ) . The fi parameter is defined as the ratio between the ^^C mass fraction and the 
corresponding solar value scaled to the initial Z Higher values of are needed to follow 
the evolution of low Z models up to the AGB tip, due to the lower initial carbon mass 
fraction and the more efficient third dredge-up. 

Our computations start from a proto-star, namely a model in hydrostatic equilibrium, 
chemically homogeneous, fully convective and cooler than the threshold for any nuclear 
burning. We follow the evolution from the pre-main sequence up to the end of the AGB 
phase, passing through the core He-flash. No core overshoot has been adopted during core 
H- burning. Semi-convection was applied during core He-burning (see also Section [5]). We 
stop the computations when, as a consequence of the mass loss, the H-rich envelope has 



^Note that fi = 1.8 corresponds to C/0~1 
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been reduced to the minimum mass for the occurrence of the TDm (see 



J3 



Straniero et al. 



2003aJ). In Table [2] we summarize the main properties of the evolutionary sequences prior 
to the thermally pulsing AGB phase. For each metallicity and initial He content we report 
from left to right: the initial mass (solar units), the core H-burning lifetime (in Myr), 
the maximum size of the convective core developed during the main sequence (in solar 
units), the He surface mass fraction after the occurrence of the First Dredge- Up (FDU), the 
luminosity at the RGB tip (logarithm of L in solar unit), the mass of the H-exhausted core 
at the beginning of core-He burning (M^g, in solar units), the core He-burning lifetime (in 
Myr), the mass of the H-exhausted core at the end of core He-burning (M^ , in solar units) 

. When comparing these models with those already 



published in 



Dommguez et al. 



(119991 ) (which where obtained with an older release of the 



FRANEC code) variations typically lower than 10% are found. These variations are due to 
the many upgrades of the input physics and to the differences in the initial composition 
(mainly the He content). For a more complete description of the pre- AGB evolution a nd 



the relevant theoretical uncertainties, the reader may refer to 



Dommguez et al. 



mm . 



2.1. The Thermally Pulsing AGB phase 

In Tables [3l HI |5] we report quantities characterizing the evolutionary sequences durin; 
TP-AGB phase for three key metalhcities: Z=0.02, Z=0.006 and Z=0.001 For each 



^The minimum envelope mass for TDU occurrence is Menv ~ (0.4 O.5)M0 at large 
metalhcities and Menv ~ (0.3 0.4)Mq at low metalhcities. 

^We define the first TP when the luminosity generated by So; reactions exceeds IQ^Lq for 
the first time. 

^A pdf file with Tables corresponding to all the computed metalhcities is available on the 
web page of the database. Note that some of the models with M = 2Mq have already been 
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model, we report (left to right): the progressive TP number (nxp)!!!, the total mass (Mxot), 



the mass of the H-exhausted core (MhJ^, the mass of the H-depleted material dredged- up by 
the TDU episode which follows a TP(AMtdu)) the maximum mass of the convective zone 
generated by the TP (AMcz), the growth of the H-exhausted core during the interpulse 
period preceding the current TP(AMh), the overlap factor the A factoil^, the duration of 
the interpulse period preceding the current TP (Atjp), the maximum temperature attained 
at the bottom of the convective zone generated by the TP (Tmax), the surface metallicity 
attained after the TDU (Zgurf), the corresponding C/0 ratio and the average number of 
neutrons captured per initial iron seed nucleus during the radiative ^^C pocket burning 
{ric). Lastly, we report the total mass of the material cumulatively dredged-up (M^^u) 
the duration of the C-rich phase of the model (tcJ^ for each model. 

The main properties of the TP-AGB phase of 2 Mq stars as a function of the metallicity 
have been extensively presented in Paper I. In this paper we explore in greater detail 
the evolution of low mass stars by discussing how the physical conditions favoring the 
AGB nucleosynthesis scenario change as a function of the initial mass. One of the most 
important quantities driving the nucleosynthesis in low mass AGB stars is the amount of 
^^C available in the pockets, which form in the transition layer between the H-rich envelope 
and the H-depleted core after each TDU. In our models, the formation of the ^^C-pocket is a 



presented in Paper I but, for a better understanding, have also been reported in this paper. 
^Negative numbers correspond to TPs not followed by TDU. 
^At the onset of the thermal pulse. 

^Namely, the overlap, in mass, between the convective zones generated by the present 
and the previous TP. 

"^"Namely, the ratio between AMtdu and AMh. 

^^That is the time spent of the AGB with a surface C/0>1. 
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X(i3c )=13*(X(i3C)/13-X(i4N)/14) 




0.001 



M/M^ 

Fig. 1. — Mass extension of the effective ^^C in the 1**, 6*^^ and 14*'* (last) pocket of the 2 Mq 
model with Z — 6 x 10"^. The pockets have been shifted in mass in order to superimpose 
their external borders. The point of the abscissa is arbitrary. 
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by-product of the introduction of an exponentially decaying profile of convective velocities at 
the base of the envelope (see Formula 1 in Paper I). This has been done in order to stabilize 
this convective boundary during TDU episodes (see Paper 1). Such a discontinuity originates 
from the opacity difference between the H-rich envelope and the He-rich intershell (i.e. the 
region between the He-shell and the H-shell, hereafter He-intershell) . Even if this additional 
mixing always works during the TP-AGB phase, the deriving opacity-induced extra-mixed 
zone is appreciable only during TDU episodes (and thus negligible during the interpulse 
phases). At the beginning of the interpulse, the few protons in this transition layer are 
captured by ^^C nuclei, producing ^^C (via the ^^C(p,7)^^N(/3"^z/)^'^C nuclear chain). A 
further proton capture may eventually lead to the production of ^^N. Consequently, after 
each TDU episode, a ^^C pocket forms, partially overlapped to a ^^N pocket. Note that 
where the abundance of ^^Ne is comparable to the ^^C abundance, a third small ^^Na 
pocket forms (see Figure 13 of Paper 1). Due to the resonant ^^N(n,p)^^C neutron capture 
cross section, ^^N is a very efficient neutron poison. However, the released protons give rise 
to a partial neutrons recycling through the following chain: ^^C(p,7)^^N(/3'^)^^C(a,n)^^0. 
Since our models are calculated with a full nuclear network, the combined effect of poisons 
and neutron recycling on the on-going s-process nucleosynthesis is properly taken into 
account. To illustrate the properties of our ^^C pockets that follow, in the following we 
use the effective ^^C fraction {i.e. the difference between the number fractions of ^^C and 
^^N in the pocket), rather than the bare ^^C amount (see the Formula in Fig. [1]). Note 
that such a quantity is related to the total number of neutrons available for the subsequent 
s-process nucleosynthesis, which takes place when the temperature in the pocket becomes 
high enough for the activation of the ^'^C(a,n)^^0 reaction. As already outlined in Paper 
I, the pockets shrink in mass ascending the AGB, following the compression of the whole 
He-intershell region. As a typical example, in Fig. [T]we plot the effective ^^C characterizing 
the I'**, e*'' and 14*^^ (last) pocket of the 2 Mq model with Z = 6 x IQ-^. The pockets 
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have been shifted in mass in order to superimpose their external borders; the point of the 
abscissa is arbitrary. The extension of the pockets decreases with the TDU number, the 
first being the largest one (AM ~9xl0~^ Mq, about 6 times greater than the last one) . 

In Fig. [2] we illustrate the variation of the total mass of effective ^^C in each pockelH^ 
as a function of the core mass along the evolutionary sequences corresponding to the 
various initial masses. We concentrate on the selected metallicities of Table [21 E] and S] 
{Z = 2 X 10"^, Z = 6 X 10~^ and Z = 1 x 10~^, respectively). Each point corresponds to a 
single ^'^C pocket and, in turn, to a single s-process nucleosynthesis episode. In this context, 
the relative contribution of a single 5-process episode to the overall s-process nucleosynthesis 
is proportional to the quantity reported on the y-axis of this plot. As extensively discussed 
in Paper I, the first few pockets dominate (qualitatively and quantitatively) the overall 
nucleosynthesis of the elements during an AGB evolutionary sequence. In practice, the 
composition in the He-intershell is freezed after just 4 or 5 thermal pulses followed by TDU 
episodes, being only marginally modified by the late ^^C burnings. In Paper I we outlined 
that, in the 2 Mq models, the first ^^C pocket is only marginally consumed before the 
onset of the following TP and thus some ^'^C is engulfed in the convective shell. We confirm 
this finding also for the other masses. This phenomenon is particularly strong at large 
metallicities {Z > 1 x 10^^), where up to 80% of ^^C in the first pocket burns convectively. 
Moreover, for the more massive and metal-rich models, a non negligible ^^C percentage 
(<30%) in the pocket formed after the 2^^*^ TDU is also ingested in the TP. Lowering the 
metallicity, the amount of ingested ^'^C progressively decreases, being completely burnt 
in radiative conditions at Z = 1 x 10^^. The convective ^^C burning has only minor 
consequences on the chemical evolution of the models (see Section [3]). Note that even for the 
first two pockets a non negligible amount of ^^C (at least 20%) burns radiatively during the 



That is the ^'^C mass integrated over the pocket region. 
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Fig. 2. — The mass of effective ^^C in the pockets formed after each TDU episode, 
SM(-'^'^Ccfr), versus the core mass. Symbols correspond to different initial stellar masses, 
as labelled in the figure. The top panel refers to models with Z — 0.02, the middle panel 
with Z = 0.006 and the lower panel with Z = 0.001. 
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interpulse phase (the exact value depending on the mass and the metalhcity). Since the first 
two pockets are the largest ones, they provide an important contribution to the s-process 
nucleosynthesis even if the amount of ^'^C burnt in radiative conditions is reduced with 
respect to the following pockets. In Fig. |2]we take into account the partial ingestion of the 
first pockets by subtracting the engulfed ^'^C in the derivation of the effective ^^C. Changing 
the mass and the metallicity of the model, a certain spread in the mass of the effective ^^C 
is found. Nonetheless, this occurrence does not necessarily imply an equivalent spread in 
the final surface s-process distributions, which basically depend on the neutrons per iron 
seed in the pocket (see Section 12]). In this connection, an interesting quantity is the average 
number of neutrons per seeds in each pocket, as reported in the last column of Tables |3l 
m [51 Recalling that the first few pockets dominate the overall s-process nucleosynthesis, 
we see that the neutrons captured per iron seed substantially increases by reducing the 
metallicity, while, for a fixed Z, it is almost insensitive to a change in the initial mass. 

Due to the partial overlap of the convective zones generated by the recurrent TPs, 
the s-processed material accumulates within the He-intershell. This occurs for the first 
convective shells, in which the production of the s-process elements within the ^'^C pockets 
compensate the dilution over the ashes of the H-shell, which are basically s-process free (see 
Fig. 16 of Paper I)Q. Afterwards, the dilution effect dominates and in the following TPs 
the heavy elements abundances in the He-intershell decrease. Fig. [3] illustrates the ratios 
between the amount of dredged-up material per TDU and the residual envelope mass as a 
function of the core mass, for the different stellar initial masses. Each point corresponds 
to a single TDU episode. Since the variation of the mass fraction of any isotope in the 
envelope as a consequence of a certain TDU episode is proportional to AMtdu/^enV) such 



^■^Actually, this material also presents an s-process enrichment, due to previous TDU, but 
it is definitely lower than the one characterizing the He-intershell. 
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Fig. 3. — As in Fig. [21 but for the ratio between the mass of H-depleted dredged- up material 
and the envefope mass. 
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a ratio well represents the efficiency of the various TDU episodes. It generally increases 
during the AGB evolution, with a sharp drop to when the TDU ceases. This ratio is also 
generally larger at low Z because the TDU is deeper. Note that the 3 Mq model with 
Z = 0.001 shows a substantially smaller AMtdu/^env than other masses with the same 
metallicity. This model develops a larger core mass and has a lower AMtdu? due to the 
shrinking of the He-intershell. Moreover, it has a larger envelope mass. Its structure is very 
similar to the one characterizing Intermediate Mass Stars (IMS) with higher metallicities. 
The modelling of IMS is made difficult from the lack of well constrained (and direct) 
observational counterp arts, even if new interesting data have been recently presented by 



Roederer et al. 



(120111 ). Moreover, the computed surface compositions of IMS strongly vary 
from author to autho r, depending on the adop ted mass-loss rate, convection efficiency and 



opacity set (see, e.g.. 



Ventura fc Marigo 



20091 ). Unfortunately, these quantities cannot be 



fixed one to the time, since they are highly degenerate. 

As already mentioned in the Introduction, the activation of a second neutron burst, due 
to the ^^Ne(a,n)^^Mg reaction, depends on the temperature attained at the base of the 
convective zone generated by a TP. This is illustrated in Fig. HJ For each initial mass a 
general trend can be seen. The maximum temperature within the convective zone generated 
by a TP progressively gro ws, reaches a ma x imum and then slightly decreases (during the 



last TPs). As discussed in 



Straniero et al. 



f l2008a 



this temperature trend depends on the 



core mass, the envelope mass and the metallicity. The maximum temperature is higher for 
larger core masses and larger envelope masses, while it is lower for higher Z. As a results, 
at Z=0.02, all stars with M < SM© barely attain the threshold of the '^'^Ne{a,ny^Mg 
reaction. On the contrary, at Z=0.001, in all stars with M > 2Mq the second neutron 
burst plays a relevant role in the resulting s-process nucleosynthesis. As we will show in the 
next section, in these models there is a substantial activation of some 5-process branchings 
and the corresponding increased production of branching- dependent isotopes (such as ^°Fe, 
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Fig. 4. — As in Fig. [21 but for the maximum temperature attained at the bottom of the 
convective zone generated by a TP. 
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^^Kr, ^^Rb and ^^Zr). In addition, light-s (Is) elemento are efficiently produced as a 
consequence of the rapid mixing taking place in the convective shell. On one hand, this 
mixing continuously provides fresh Fe seeds in the innermost (and hotter layers) of the 
convective zone. On the other hand, it prevents the pile up of the synthesized material, 
which is required to build up the heaviest species (as it occurs during the radiative ^^C 
burning within the pockets). 



3. Nucleosynthesis 

Three main peaks are expected in the s-process distribution. They correspond to nuclei 
with a magic number of neutrons (50, 82 and 126), whose binding energy is particularly high 
and, in turn, the neutron capture cross-section is particularly low. The average abundances 
of s-process elements around N=50 and N=82 represent respectively the so-called Is (light 
s) and hs (heavy s) indexes. As in Paper 1, we define: 

• [Is/Fe] = ( [Sr/Fe] + [Y/Fe] + [Zr/Fe] ) /3; 

• [hs/Fe] = ( [Ba/Fe] + [La/Fe] + [Nd/Fe] + [Sm/Fe] ) /4. 

The third peak (at N=126) is represented by lead ([Pb/Fe]). Combining these quantities, 
the spectro scopic inde x es [hs /ls] , [Pb/ls] and [Pb/hs] can be derived!^. Basing on simple 



arguments. 



Cameron 



( 119571 ) demonstrated that these spectroscopic indexes basically 
depend on the number of neutrons available per iron seed for the first time. In stellar 
interiors, however, the physical and chemical conditions may be more complex. In the case of 



^^Those belonging to the first peak of the s-process (Is, light elements: Sr-Y-Zr). See next 
Section. 



^^We adopt the usual spectroscopic notation: [a/b]=log(a/b)-log(a/b) 
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the radiative ^^C burning taking place during the interpulse period in low mass AGB stars, 

the number of neutrons corresponds to the number of ^'^C nuclei in the pocket. However, 

most of these neutrons are captured by light elements (poisons) and do not contribute to 

the s-process. As already specified in Section I2.H ^"^N (the most abundant poison in our 

models) partially recycles neutrons. Nevertheless, as a first order approximation, we may 

assume that the number of neutrons per seed available for the s-process is given by: 

n{neutrons) — n{poisons) n{^^C) — n("'^^N) 
n{seeds) n(^^Fe) 



The numerator of this equation coincides with the effective ^^C introduced in Section flA\ 
Within the ^'^C pocket it varies according to the protons profile left by the TDU. Note that 
both the ^^C and the ^^N in the pocket share the same (primary) origin, since both are built 
from the primary ^^Ccj. On the contrary, the iron abundance within the pocket depends on 
the initial stellar metallicity. 

The behavior of the three s-process peaks at the last TDU as a function of the 
metallicity and for different masses is illustrated in Fig. [5l For selected masses and 
metallicities, the corresponding values are reported in Table |6l More complete Tables can 
be downloaded from the FRUITY web-based platform. We find that [Is/Fe] is barely flat 
for all the metallicities, the only 1.5 Mq exhibiting a mild increase when decreasing the 
metal content. As to [hs/Fe], all masses show a steep increase with the metallicity. This 
behavior is even more evident for [Pb/Fe], which peaks at the lowest computed metallicity. 
Given the steep increase of the neutrons per iron seed with decreasing metallicity (see 
column 13 in Tables [3l |l]and|5]), the heaviest species are easier synthesized at low Z. In 
Fig. [6] we report the s-process spectroscopic indexes [hs/ls] (upper panel) and [Pb/hs] 
(lower panel). The first quantity continuously grows with decreasing metallicity, reaching 



Primary elements are direct by-products of H or He burning. 
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Fig. 5. — Behaviors of the three s-process peaks ([Is/Fe], top panel; [hs/Fe], middle panel; 
[Pb/Fe], lower panel) as a function of the metallicity predicted by our models for different 
initial masses. 
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Fig. 6. — As in Fig. [5l but for the s-process indexes [hs/ls] (upper panel) and [Pb/hs] (lower 
panel) . 
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a maximum at Z ~ 3 x 10~^ (eventually decreasing at lower metallicities), while [Pb/hs] 
shows a flat profile from Z ~ 2 x 10^^ down to Z ~ 6 x 10^'^ and then grows at the lowest 
metallicities. A full discussion on the behavior of th ese quantities as a function of the 
metal content can be found in 



Bisterzo et al. 



(120101) • As already outlined, these quantities 



mainly depend on the neutrons per iron seed and on the abundances of neutron poisons. In 
our models, the amount of ^^C in the pockets basically depends on the core mass (see Fig. 
[1]). At fixed metallicity the core masses for different models do not vary much (apart from 
the Z = 1 X 10~^ models, but see below) and, therefore, there is no particular trend with 
the initial mass. Our spread at fixed metallicity refiects the slight variation of the total 
number of neutrons per iron seed in each pocket for different initial masses. This limited 
spread derives from the fact that, in our models, only convection is considered as a mixing 
mechanism: the int rodu ction of additional m ixing phenomena (such as rotation; see e.g. 



Langer et al. 



3 



19991 and 



Herwig et al. 



20031 ). able to dilute the available ^^C on a larger 



mass region, could lead to a lower neutron to seed ratio (see Section [6^ . 

Other interesting nucleosynthetic features emerge from the analysis of light elements 
(see Table [6], in which the final C/0, C/N and ^^C/^^C ratios are also given). Among them, 
carbon and fiuorine show a clear anti-correlation with the initial metallicity. Carbon is of 
primary origin, because it is produced during TPs via the 3a reaction. The production 
of fiuorine basically depends on the availability of ^^C in the He-intershell (see Paper 
I for a detailed description of the nuclear chain leading to ^^F nucleosynthesis; see al so 



Jorissen et al. 



Abia et al. 



1992 



2010 . 



Forestini et al. 



1992 



Mowlavi et al. 



1996 



Lugaro et al. 



200 



20111 ). A non negligible amount of ^'^C is available within the ^^C-pocket 
and in the ashes of the H-burning shell. As already discussed in Paper I, the former is 
weakly dependent on metallicity. The latter directly scales with the CNO abundances in 
the envelope, which, in turn, depend on the (primary) ^^C dredged-up during all previous 
TDU episodes. The resulting fiuorine can be roughly considered a primary element. 
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The nitrogen surface enrichment is mainly a consequence of the FDU and, at the lowest 
metallicities, of the H-shell erosion by the envelope during TDU episodes. This latter 
contribution is not appreciable at higher metallicities because of the already large ^^N 
surface abundance. 

The neon surface enrichment derives by the dredge-up of ^^Ne, which is synthesized 
during TPs via the nuclear chain ^^N(a,7)^^F(/3+z/)-'^^0(a,7)^^Ne. Direct evidence of 



neon production has recently been detected in Galactic planetary nebulae (IMilingo et al. 



2OIOI ). Since the nitrogen abundance in the He-intershell depends on the quantity of 
carbon previously dredged-up (due to the conversion of ^^C to ^^N via the CNO cycle), 
^^Ne is a primary element (thus its surface overabundance is expected to increase at low 
metallicities). 

In ou r models, sodium can be synthesized via proton or neutron captures (see also 



Mowlavi 



19991 ). A sodium pocket forms via the Ne(p,7) Na reaction in the upper part of the 



He-intershell, where protons are left by the receding convective envelope after TDU episodes 
(see Paper I). Moreover, ^^Na is synthesized via the ^^Ne(n,7)^^Ne(/3~z/)^^Na nuclear chain, 
during both the radiative burning of the ^^C-pocket and the convective ^^Ne burning in the 
TPs. 

Magnesium isotopes follow different nucleosynthetic paths. ^''Mg, which is the most 
abundant magnesium isotope in the Sun, is marginally produced via neutron captures on 
^^Na during the radiative burning of the ^'^C pocket and it is weakly destroyed during 
TPs. On the contrary, ^^Mg and ^^Mg are mainly created during TPs via neutron captures 
processes and via the ^^Ne(Q;,n)^^Mg and the ^^Ne(a,7)^^Mg reactions (which contribute 
to the ^^Mg and ^^Mg nucleosynthesis, respectively). In our models we generally find a 
low final surface magnesium enrichment, with an increasing trend toward low metallicities. 
This derives from the larger relative contribution of ^^Mg and ^^Mg, due to the increased 
efficiency in the activation of the ^^Ne(Q;,n)^^Mg and ^^Ne(a,7)^^Mg reactions (see Section 
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All our models become C-rich [i.e. their final C/0>1), but the 1.5 Mq model at 
Z = 2 X 10^^. As already mentioned, non convective mixing episodes (the so-called 
extra-mixing) for both the RGB and the AGB phases are not considered in these models. 
The origin of this extra-mixing is far from being established on theoretical grounds. The 
precise estimate of its efficiency during both the giant branches can be obtained only by 



comparing parameter i zed models with obseryational constraints ( Charbonnel fc Zahn 



Eggleton et a 



Denissenkov 



2008: 



20101 



Nordhaus et a 



Palmerini et al. 



2008 



2011 



Denissenkov et a 



Denissenkov et al. 



20091 



Busso et al. 



20071 : 



2010 



20111 ). The extra- mixing is 
expected to af fect the CNO abundarices, with a substantia l reduction of the ^^C/^^C isotopic 



ratio (see also 



Lebzelter et al. 



2008 



Lederer et al. 



20091 ) and, eventually, a lowering of the 



C/0 and C/N ratios. For this reason, special attention has to be paid when comparing the 
surface abundances of these elements (or isotopes) with their observed counterparts. 

A selection of net stellar yields for different masses and metallicities is reported in 
Tables [71 [S] and [3 (full Tables are available on line at FRUITY web pages). These data 
represent the first homogeneous set of light and heavy elements theoretical yields for AGB 
stars. The completeness of our database could help to r educe the uncertainties affecting 



Romano et al. 



f!2010h . At high 



Galactic chemical evolutionary models, as underlined by 
metallicities (Z ~ 2 x 10^^), the largest absolute yields (for both light and heavy elements) 
come from the 2.5 Mq and the 3.0 Mq. At lower metallicities (Z ~ 6 x 10^'^) the principal 
contribution derives, for all chemical species, from the 2.5 Mq model. At Z ~ 10^'^, the 
initial mass of the principal polluter further decreases, presenting the largest yields in the 



2.0 Mq mass. Weighing the yie. 



one proposed by ISalpeter 
metallicities. 



ds with an Initial Mass Function (IMF), for example the 



( 119551 ). the contribution from the 2.0 Mq stars dominates at all 
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The amount of processed material from various masses at different metallicities results 
from a combination of many factors, as the core mass, the efficiency of TDU, the envelope 
mass, the number of TPs, the thickness of ^^C pockets and the surface temperature (which 
drives the mass-loss rate). As a rule of thumb, a progressive decrease of the metallicity 
leads to a reduction of the yields from stars with M ~ (2.5 ^ 3.O)M0, making stars with 
M ~ (1.5 2.O)M0 the principal polluters of the ISM. Note that, as already pointed out, 
the explicit treatment of physical phenomena currently not accounted for (such as rotation) 
could significantly modify these quantities and, in particular, the surface distributions of 
the heavy elements. 

In order to better interpret the nucleosynthetic paths described above, we highlight 
some general rules characterizing our models: 

• the main neutron source is represented by the radiative burning of ^^C within the 

pockets; 

• the nucleosynthesis occurring in the first 4—5 pockets mainly determines the final 
surface distributions; 

• the contribution of the additional ^^Ne(a,n)^^Mg neutron burst is generally negligible; 

• the efficiency of the TDU increases with decreasing the metallicity; 

• at fixed metallicity, the surface chemical enrichment basically depends on the total 
amount of dredged- up material and on the mass extension of the convective envelop J^. 

However, there are a few exceptions. As already described in Paper I, the first ^'^C-pocket of 
the 2 Mq models does not completely burn during the radiative interpulse period and some 



^ ''See Section H] for a discussion about the effects induced by different mass-loss rates on 
the final elemental surface abundances. 



-27- 



^^C is ingested and burnt convectively in the following TP. We confirm this finding for the 
other masses. This episode has no effects on the physical evolution of the TP, while it has 
interesting effects from an isotopic point of view. In particular, the large neutron densities 
attained (see Paper I) allow the nucleosynthesis of some neutron-rich isotopes, such as the 
short lived ^°Fe. 

A characteristic feature of our 3 Mp^ model at Z = 10~^ is represented by the occurrence o f 



the so-called Ho t-TDU (iGoriely fc Siess 



Lau et al. 



2004 : 



Herwig 



2004 



Campbell &: Lattanzio 



2008 



20091 ). In fact, for this model proton captures occur at the base of the convective 



envelope during the TDU phase. In Paper I, we already found the Hot-TDU to occur in the 
M=2 Mq model with Z = 10~^. As a first consequence, the model shows a definitely lower 
i2q^i3q igotopic ratio with respect to other masses (see Table |6]). Note that the occurrence 
of the Hot-TDU depends on the mixing algorithm we applied at the inner border of the 
convective envelope (see Section |2]). The calibration of its /3 parameter has been done in 
order to maximize the effective ^^C in the pockets and, consequently, the resulting s-process 
nucleosynthesis. Actually, it is not clear if the same calibration can be also applied to AGB 
stars developing a larger core mass, as our 3 Mq model at Z = 10~^. With respect to the 
other calculated masses, this model also reaches higher temperatures during TPs, which 
imply a more efficient activation of the ^^Ne(a,n)^^Mg neutron source. This leads to a non 
negligible contribution to [Is/Fe]. Consequently, the [hs/ls] index decreases, leading to a 
larger spread of this quantity when compared to the other metallicities. 
It is worth to mention that the nucleosynthesis of some isotopes (the ones whose production 
requires a high neutron density) is favored by the activation of the ^^Ne neutron source. 
In Table [10] we list some isotopic ratios involved in branchings along the s-process path. 
The dependence on the metallicity has already been discussed in Paper I, demonstrating 
that these ratios increase with decreasing metallicity, due to the increased contribution to 
the 5-process nucleosynthesis of the ^^Ne(a,n)^^Mg. The dependence of the corresponding 
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isotopic ratios on the initial mass is clearly shown in Table [101 With the exception of the 
most metal-rich models, the 3 Mq models show the largest ratios, confirming again that the 
^^Ne(Q;,n)^^Mg neutron source is particularly important for the more massive AGB stars. 



4. Theoretical uncertainties 

When dealing with observations, it is mandatory to perform a detailed discussion of 
the errors introduced by the instruments and the techniques adopted to analyze data. A 
similar procedure should be followed when computing theoretical stellar models. In this 
case, error-bars cannot be derived (due to the high degree of degeneracy among the input 
parameters), but the sensitivity of the model with respect to changes in the input physical 
parameters can be evaluated. A similar analysis was already performed in Paper I for 
a 2 Mq AGB model with Z = 10^'^. Considering the larger metallicities of the models 
presented in this work, we focus on a M = 2Mq with Z = 6 x 10~^ (hereafter FRANEC 
Standard Model, FSM). We assume as our FSM the model computed with the following set 
of parameters: 

• the free parameter of the mixing length theory is set to q;fsm=2.15 (see Paper I); 



the AGB mass-loss rate (M-^^f) is empirical 



luminosity relation (see 



Straniero et al. 



2006 



y calibrated on the observed period- 
for details); 



• the free parameter that determines the decreasing profile of the convective velocities 
at the base of the envelope is set to /3fsm=0.1 (see Paper I). 



In order to estimate the sensitivity of our models to the afore-listed inputs we calculate six 
test models. We vary the mixing length parameter a, the mass loss law and the /3 parameter 
which governs the efficiency of the velocity profile at the border of the convective envelope: 
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1. ai = apsM + 0.35 = 2.5; 



2. a2 = apsM — 0.35 = 1. 



3. M,^«^ = M^g! * 2; 



5. A = /3fsm + 0.02 = 0.12; 



6. 132 = /3fsm - 0.04 = 0.06. 



A variation of the mixing length parameter a affects not only the temperature run in the 
more external layers of the star, where the local temperature gradients largely deviate from 
the adiabatic value, but also the velocity prof ile of the convective bubble in the whole 



convective envelope (see 



Piersanti et al. 



20071). 



We evaluate the effects of a variation of the AGB mass-loss rate by multiplying or dividing 
its efficiency by a factor 2. At fixed pulsational period, the resulting theoretical mass loss 
laws reproduce the m inimum and the maxim um mass- loss rates observed in Galactic AGB 



stars (see Figure 3 of lStraniero et al. 1120061 ). 

Finally, we present two models with a different /3 parameter in order to estimate the effects 
induced by a change in the slope of the exponentially decreasing profile of convective 
velocities at the base of the envelope. We choose values leading to a reduction of the 
effective ^^C within the pockets by roughly a factor 4 with respect to our FSM (see Figure 
7 of Paper I). 

In Table [TT] we report, for our FSM, the final core mass (Mh), the mean bolometric 



magnitude of the AGB C-rich phase (M^^jJ^, the duration of the AGB C-rich phase (r^) 



^'^Calculated as the average luminosity weighted over the time spent by the star on the 
AGB when C/0>1. 
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the ^^C yield (as representative of primary-like light isotopes), the ^^Y, ^^^La and ^osp]^ 
yields (as representative of the three s-process peaks) and the final [hs/ls] ratio. For the 
test cases, we list normalized percentage differences with respect to the FSM. 
The physical quantities listed in Table [11] (final core mass, mean bolometric magnitude of 
the C-rich phase and duration of the C-rich phase) mildly depend on the choice of the input 
parameters. Maximum variations are below 30%. This makes our theoretical models quite 
robust when comparing with the observed counterparts (see Section [6?T|) . 
The test with an enhanced a parameter and the test with a reduced mass loss rate present 
similar trends: the mean bolometric magnitude of the C-rich phase decreases (and, thus, 
the mean luminosity increases), while the C-rich phase duration and the chemical yields 
increase. The reason is twofold. In fact, a higher a parameter implies more efficient TDUs 
and a slightly higher surface temperature, which translates into a milder mass loss rate. An 
opposite trend is found for the test with a reduced a parameter which, thus, resembles the 
test with an increased mass loss rate. Note that, for the aforementioned test models, the 
heavy elements yields roughly scale as the light ones and the final chemical patterns show 
minor changes with respect to the FSM (see their [hs/ls] ratios). 

A different prescription for the mass loss law leads to negligible differences in the s-process 
abundances because the first ^^C pockets largely determine the final heavy elements 
enrichment in the envelope. For the same reason, minor changes in the relative abundances 
of s-process elements are found in the a-varied tests, even if a variation of this parameter 
leads to different TDU efficiencies. However, a further aspect has to be analyzed. Within 
the framework of the mixing length theory, the value of the convective velocities directly 
scales with a. Thus, a variation of such a parameter should also modify the efficiency of 
the mixing algorithm we adopt at the base of the convective envelope. However, we verified 
that this has minor effects on the formation of the ^^C pockets, with variations lower than 
15% in the effective ^^C. This is not the case for the two tests with a different /3 parameter. 
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In fact, these models present complete different heavy elements distributions with respect 
to the FSM. In both cases there is a sensible decrease of the absolute surface abundances 
(as proved by the yields) and a change in the final shape of the s-process distribution (as 
showed by the [hs/ls] ratios). This is due to the fact that the amount of neutrons available 
for the nucleosynthesis of heavy species within the ^^C pockets decreases in both cases. 
The reasons are, however, different. In the /3=0.06 case, the total amount of ^^C within 
the pockets strongly decreases due to the steeper decreasing profile of convective velocities. 
Instead, in the (3=0.12 case the ^^C pocket is completely overlapped with the ^^N pocket. 
This acts as a neutron poison via the resonant ^^N(n,p)^'^C reaction, thus penalizing the 
heavy elements nucleosynthesis (see Section [3]). 



5. Comparison with other stellar models 



In previous Sections we discussed in detail our theoretical models, by analyzing their 
physical and chemical properties as a function of the initial mass and metal content. 
Moreover, we evaluated the uncertainties deriving from different choices of the input 
parameters. In this Section we present a comparison between our results and other 
evolutionary stellar models available in the literature. A similar comparison is in general 
not straightforward, since the various codes significantly differ in many of the adopted 
input physics (treatment of convection, definition of convective borders, mass-loss law, 
initial chemical abundances, opacities, equation of state and nuclear network). As already 
discussed in Section El during the pre -AGB phase our models are very similar to those 



presented in 



Dommguez et al. 



different codes. As discussed in 



9991) . We refer to that paper for a comparison between 



Straniero et al. 



( l2003al ). a key physical quantity to 



understand AGB evolution is the core mass at the beginning of the TP-AGB phase (we 
identify it as the core mass at the first TP in which the helium luminosity exceeds lO^L©, 
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TP1^ 



In Fig. [7] we plot our M^^^ for three select ed metal. 



and Z = 0.003) arid corn pare them to models from 



Karakas 



icities (Z = 0.02, Z = 0.008 



Weiss &: Ferguson 



20091) 



(120031 ) (hereafter K03) and 



(hereaft e r WF 09). For the M = ?, Ma, model with Z = 0.0 2 



Herwig 



fl2000l ) (hereafter HeOO) and 



Stancliffe et al. 



(120041 ) 



we also report data from 
(hereafter Sta04). 

The value of M^^^ basically depends on the core mass attained at helium ignition (M^^ 
in Table [2]) and on the duration of core He-burning. M^g is built during the RGB phase 
for an initial mass lower than 2 Mq. M^g depends on several input physics, like the 
neutrino cooling rate, the ^^N(p,7)^^0 reaction rate and the equation of state for the 
degenerate matter of the RGB core (including the Coulomb correction). For higher initial 
masses, which develop a sizeable convective core during central H-burning, the treatment 
of convection (overshoot or not) may also affect M^^. On the other hand, the duration 
of the core He-burning phase depends on the treatment of convection (semi-convection 
or mechanical overshoot) and on the adopted values for the key reaction rates 3a; and 
^^C(a,7)^®0. In addition, we stress that the difference in the initial H and He abundance 
has a significative influence on the whole stellar model evolution. 

Our models are in fair agreement with those presented by K03, who operated similar choices 
(no core overshoot during core H-burning and semi-convection during core He-burning). 
A similar M^^^ is also attained by Sta04, who also did not apply convective overshoot. 
Note that Sta04 simultaneously solve the equations of the chemical evolution (including the 
convective mixing) and those of the physical evolution. Such numerical algorithm could 
mimic a treatment of semi-convection similar to the one we adopt. The reason of the 
notable difference between our models and the WF09 ones for the M = 2Mq masses is not 
clear. We note that, for these masses, the duration of their core He-burning phase is larger 
than ours. Such an occurrence suggests that their M^g should be lower than ours (lower He 
core masses imply lower luminosities and then larger He- burning lifetimes). Unfortunately, 
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WF09 did not report M^^. On the contrary, their M = SMq models have larger M^^^ than 
us and K03. This could derive from the application of convective overshoot during central 
H-burning. Note that an even larger M^^^ is attained by HeOO, who also applied overshoot 
during both core hydrogen and helium burnings. 

In order to compare the TDU efficiency among different codes, in Fig. [8] we report the 
variation of the A factor as a function of the core mass for our AGB model with M = 3Mq 



and Z=0.02. We pl ot the same quantity 



the FRANEC code (IStraniero et al. 



authors (HeOO; Sta04; 



br a model computed with an older version of 



19971) (herea fter Stra97) and for calculations by other 



Karakas fc Lattanzio 



20071 . hereafter K07). The TP-AGB phase of 



the Stra97 model was computed without the exponentially decaying profile of convective 
velocities at the base of the convective envelope (see Section [2]). As a by-product, the TDU 
efficiency increases (in fact the maximum A grows from their 0.4 to the current 0.6) and, 
in addition, the TDU occurs earlier (see also the K07 model ). Further differences ca n be 



Straniero et al. 



fl2003ah . the 



ascribed to the adopted mass-loss law. In fact, as shown by 
TDU efficiency, and then A, substantially depends on the residual envelope mass. The K07 
model is calculated with a milder mass-loss rate than the one adopted by Stra97 (Reimers' 
formula with ri=3 during the AGB phase). Thus, for the same core mass, their model has a 
larger envelope mass and, in turn, larger A values. For the same reason, the Sta04 model, 
which is calculated without mass-loss, attains larger A values than ours. Note that the 
most efficient TDU is obtained by HeOO, even if he included a moderate mass- loss (Reimers' 
formula with ?7=1 during the AGB phase). However, in his model the convective overshoot 
is also applied at the base of the conve ctive shells, this fac t causing stronger TPs and. 



consequently, more efficient TDUs (see 



Lugaro et al. 



20031 ) . Fin ally, further differences 



could be ascribed to th e choice of the ^'^N(p,7)^^0 reaction rate ( IStraniero et al. 



2000 



Herwie &; Austin 


2004 




The 


Forestini & Charbonnel 



(119971) and the 



Karakas 



feoioh (hereafter KIO) are the 
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Fig. 8. — A factor as a function of the core mass for an AGB model with M = SMp) and 
Z=0 .02 as reported by d ifferent authors (FRUITY: th is pa per: Stra97: IStraniero et al. 
1997 : HeOO: Herwig 2000 : Sta04: Stanchffe et al. 20oi K07: Karakas fc Lattanzio 2007 ). 
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only papers we found in the literature in which yields from full stellar evolutionary 
AGB models have been discuss ed ( we exclude sy nthetic yields from our analysis; see e.g. 



Groenewegen &: De Jong 



1994 and 



Marigo 



200 ll ). We compare the yields of our M = 3M( 







and Z=0.02 model with the one presented by KIO, by considering her model including a 
partial mixing zone of 2 x 1O~^M0 after any TDU episode. For its physical characteristics, 
it is the most similar model to ours; as a matter of fact, the cumulatively dredged- up 
masses are nearly the same (Mf^^^ = 7.93 x IQ-^Mq; M^g^j^^^ = 7.54 x lO-^M©, see 
Table [3]). Note that while we have computed a full network stellar model calculations, the 
yields presented by KIO are obtained by means of a post-process calculation. In Table [12] 
we report a selection of net stellar yields. The analysis is limited to light isotopes since 
heavy elements are not included in KIO calculations. For each isotope we report our net 
yield (labelled as FRUITY), the corresponding value tabulated in KIO and the percentage 
difference. In spite of the many physical differences, the resulting yields show a reasonable 
agreement. The differences concerning ^*^Ne, ^^Mg and ^®Si basically derive from the small 
AGB net production (or destruction) as compared to their initial abundances. 



6. Comparison with observations 

In the previous Section, we presented a comparison among our models and other 
models available in the literature. As we have shown, the different results can be ascribed 
to the different adopted input physics. This comparison, however, does not indicate the 
most reliable input physics or numerical scheme to be preferred. In order to do that, we 
have to be guided by observations. This is the purpose of the present Section, where we 
compare our theoretical models to observational counterparts. Here we concentrate on the 
Luminosity Function of Carbon Stars and the s-process spectroscopic indexes at different 
metallicities. 
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6.1. Luminosity Function of Carbon Stars 



As already outlined in previous Sections, the majority of our TP-AGB models present 
final C/0 ratios larger than 1; thus, their observational counterparts are C-rich stars. The 
Luminosity Function of Carbon Stars (hereafter LFCS) represents a good test indicator for 
our theoretical prescriptions, given that it links a physical quantity (the luminosity) with 
the chemistry of the model (its surface carbon abundance). By comparing our theoretical 
LFCS with observational data, we can simultaneously check if the final structure (basically 
the core mass) and the surface chemical abundances of our models (which depend on the 
TDU efficiency and the mass loss law) are reliable or not. 



Iben 



( 1l98ll ) firstly highlighted the incapability of theoretical AGE models to reproduce the 
LFCS of Magellanic Clouds. In his models, TDU was found to occur for core masses larger 
than 0.6 Mq (thus corresponding to AGE stars of Intermediate Mass: 4 < M/Mq < 7). 
However, in these objects the temperature at the base of the convective envelope is high 



( Sueimoto 


1971; 


Iben 


1973; 


Karakas 


2OO3I: 



Ventura fc D'Antona 



20051 ). basically 



prevents these stars from becoming C-rich, due to the conversion of the dredged-up carbon 



into nitrogero. Therefo re, a discrep ancy be tween theo r y and obse rvations arose, c a. 



"carbon stars mystery" (jiben 



198l|). Later, 



Lattanzio 



fll989h and 



Straniero et al. 



led 



(|l997|) 



found TDU to occur in stars with masses as low as 1.5 Mc^, thus allowing th e possibility 



Izzard et al. 



mm concluded 



to form C stars at lower luminosities. Notwithstanding, 
that full stellar evolutionary models cannot reproduce the Magellanic Clouds LFCS, being 
the evolutionary models too luminous with respect to the observed distributions. These 
authors concluded that TD U must operate earlier on the AGE or at lowe r core masses. 



More recently, two groups (jStanclifi^e et al. 



2005 



Marigo fc Girardi 



20071 ) stated to have 



19 



Note that none of the models presented in this paper suffers for Hot Eottom Eurning. 
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Stancliffe et al. 



(|2005[ ) achieved a good 



reproduced the LFCS of Magellanic Clouds 
agreement for the Large Magellanic Cloud and suggested a shift toward lower masses for 
stars experiencing TDU in order to reproduce the LFCS of the Small Magellanic Cloud. 



A ful 



Marigo fc Girardi 



match for both the Magellanic Clouds LFCS has been found by 
( l2007l ). However, these authors used synthetic stellar models and their agreement derives 
from the ad-hoc calibration of stellar properties during the TP-AGB phase (in particular 
the onset and efficiency of the TDU process}^ 



A different approach has been followed by 



Guandalini et al. 



(|2006[ ). who performed an 



analysis of Galactic carbon stars at infrared (IR) wavelengths. They suggested that the 
inconsistency between theory and observations is due to the underestimation of the IR 
contribution of these objects. In fact, since these stars are very cool, they emit a large part 
of their radiation at IR wavelengths. Actually, these authors found a distribution peaked at 
Mboi = —5, thus in better agreement with model predictions, with a non negligible luminous 
tail up to Mboi = —6. It has to be stressed that their sample contains a large fraction of 
stars (about 45%, see their Table 3) whose distance is poorly known and should be better 
determined. The use of recently improved observational data (in particular distances) could 
lead to significant changes in their observed LFCS. 

In Fig. |9] we report our theoretical L FCS (solid histogram) a nd a revised version of the 



observed Galactic LFCS presented by 



Guandalini et al. 



2006 



(dotted histogram; courtesy 



of R. Guandalini). The y-axis (fractional numbers) is defined as the fraction of carbon 
stars in the bin of luminosity with respect to the total number of carbon stars. The 
observed LFCS consists of 111 objects, whose estimate of absolute luminosity can be 
obtained in two ways. Fir stly, thanks to the use of Hipparcos parallaxes recently revised 



(Ivan Leeuwen et al. 



20071 ) and of IR data from the ISO and MSX missi ons or, secondly, 



from the pe riod-Luminosity re . 



(120061) (see 



Guandalini et al. 



atioii for Mira C-rich sources publish ed in 



Whitelock et al. 



2006 



and 



Guandalini fc Busso 



2008 



for the technical 
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Fig. 9. — Theoretical (solid histogram) Luminosity Functio n of Galactic Carbon Sta rs com- 
pared with an upgrade of the observed one presented in ICuandalini et al. I l2006l (dotted 
histogram; courtesy of R. Guandalini). 
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procedure followed to construct the LFCS). The theoretical LFCS has been obtained by 
assuming that the current AGB population of the Galactic disk consists of stars born 
at different ages (and thus with different masses) and with different metallicities. We 
hypothesize that the disk formed 10 Gyr ago and, increasing the age with time steps of 
10^ yr, we evaluate the number of AGB stars currently laying on their TP-AGB during 
each time step. We calculate the lower and upper mass on the AGB by interpolating the 
physical inputs (such as the main sequence lifetime, the AGB lifetime and the bolometric 
magnitudes along the AGB) on the grid of computed models. The upper mass has been 
fixed to 3.0 Mq, while the minimum mass depends on the metallicity. In order to cover 
the low mass tail of the LFCS, we calculate additional models with a restricted nuclear 
network. At high metallicity {Z = 2 x 10~^ and Z = Zq) the minimum initial mass for a 
star to become carbon rich is 1.5 Mq, then decreases to 1.4 Mq aX Z = 10^^, to 1.3 Mq at 
Z = Qx 10"^ to 1.2 Mq at Z = 3 X 10"^ and,^ finally, to 1.1 Mq at Z = 10_VWe adopt a 



linear decreasing Star Formation Rate (SFR)|^ and the IMF proposed by 



Salpeter 



fll955h . 



Moreover, we assume a theoretical flat metallicity distribution (see below). 
Our curve matches the observed one, show ing, within observation al errors, the same 



peak around Mboi ~ —5.0. Note that, unhke 



Stancliffe et al. 



(120051 ) ■ we do not apply any 



normalization coefficient to match the peak of the corresponding observations. Moreover, 
we highlight that the choice of a different IMF, SFR or metallicity distribution, does lead to 
very similar results. In the upper panel of FigJTOlwe report our standard luminosity function 
(solid curve), a test case with a Salpeter's IMF law with a=1.35, a test case computed with 
a constant SFR over the whole Galaxy evolution (dotted curve) and a test with a metallicity 
distribution as derived from observations of local G-dwarfs (Dual Metallicity Distribution, 



^''We assume an initial amount of gas converted to stars of 8 Mq per yr and then we apply 
a linear decrease down to a current value of 3 Mq per yr. 
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Fig. 10. — Upper Panel: theoretical Galactic LFCS calculated with different input param- 
eters: Standard Case (solid curve), modified Salpeter's IMF (short-dashed curve), constant 
SFR (dotted curve) and metaUicity distribution as observed in local G-Dwarfs (long-dashed 
curve). Lower Panel: LFCS for selected key metaUicities: Z — 1 x 10~^ (dotted curve), 
Z — 6 X 10~^ (dashed curve) and Z = Zq (solid curve). 



-42 - 



DMD; iRocha-Pinto &: Maciel Ill996l ). As already stated, the differences between different 
cases are negligible. Moreover, we evaluate the weight that each metallicity has on the 
final luminosity function. In the lower panel of Fig. we plot some LFCS corresponding 
to selected key metallicities: Z = 1 x 10~^ (dotted curve), Z = 6 x 10^^ (dashed curve) 
and Z = Zq (solid curve). It emerges that larger is the metallicity, thinner is the LFCS, 
being the Mboi peak always around —5.0 . The flattest LFCS is obtained at Z = 1 x 10^^, 
due to the decrease of the stellar mass for the occurrence of TDU (which accounts for the 
faintest stars) and to the large increase of the final core masses for stars with initial masses 
M > 2.5Mq (which accounts for the brightest stars). 

More stringent constraints would come from the comparison of theoretical LFCS with 
the ones derived in the Magellanic Clouds (MCs). This basically derives from a more precise 
determination of MCs distances with respect to Galactic carbon stars. Notwithstanding, we 
must remind that the photometric measurements of C-rich stars belonging to MCs are often 
incomplete, since they lack of the contribution from the thermal IR band (A > 2/im). As a 
matter of fact, these objects emit the majority of their flux (~90%) at these wavelengths, 
while the most widely adop ted data used to built up the LFCS were bas e d on optical 



(v isible or near-IR data) 



or 



Costa fc Frogel 



Groenewegen et al. 



by 



Westerlund et al. 



1996h for L MC. Recently, 



fl2007h and 



1986) and 



Rebeirot et al. 



Ziilstra et al. 



Lagadec et al. 



(120061) 



1993) for SMC, 



Sloan et al. 



fl2006h 



(120071 ) report the luminosities of a sample 



of about 60 C-stars belonging to MCs, as derived by combining ground-based observations 
in the near-IR (J,H,K and L) with Spitzer Space Telescope spectrosco pic observations. In 



the up per panel of Fig. [TT] we compare th e LMC LFCS as reported by 



(120071 ) (dotted curve) and IZijlstra et al. I (120061 ) (dashed curve). Although the two 



Groenewegen et al. 



observational works share the same data set, the resulting Luminosity Functions appear 
rather different. Nevertheless, they both peak at Mboi = —5, thus in good agreement with 
our theoretical LFCS. This latter curve have been constructed by assuming a constant SFR 
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Fig. 11. — Upper Panel: t heoreti cal and obseryationa l LFCS for the LMC system (Groe07 



Groenewegen et al. 



2007; Zi06: 



Ziilstra et al. 


2006 


); I 


jower 


em (S106: 


Sloan et al. 


2006 



lOwer Panel: theoretical anc 



Lagadec et al. 



obser- 



20071) 



See text for details. 
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and 3 equally weighted metallicities (Z = 3 x 10 ^, Z = 6 x 10 ^ and Z = 8 x 10" 



observed dis tributions ar e 
sample from ISloan et al. 



Tom 


Groeneweeen et al. 


(2006 


) and 


Lapadec et al. 



fl2007h 



(dotted curve) and a combined 



(12007h (dashed curve). For SMC, the 



observational LFCS shows a larger spread, but the average Mboi is similar to the LMC one. 
The lower SMC mean metallicity forces us to use models with different Z to construct our 
theoretical LFCS (Z = 1 x 10-^ Z = 3 x 10"=^ and Z = 6 x lO'^). As already explained, 
this makes our distribution flatter and, thus, in quite good agreement with observations. 
The statistical samples used for the observed CSLFs reported in Fig. [TT]are quite poor, but 
we are confident that future developments of ground and space-based IR-facilities would 
greatly improve the observational constraints to AGB theoretical models. 



6.2. s-process indexes 

In this section we compare the predicted spectroscopic indexes [hs/ls] and [Pb/hs] 
to the observed composition of various s-enhanced stars. Among these stars, we should 
distinguish between AGB stars [i.e. those presently undergoing thermal pulses and TDU 
episodes), post-AGB stars and less-evolved stars (dwarfs or giants) belonging to binary 
systems. In the latter case, the s-enhancement is likely due to the pollution caused by the 
intense wind of an AGB companion (hereinafter extrinsic s-enhanced stars). 
In principle, such a comparison requires the knowledge of the evolutionary status of the 
star, the initial mass and the metallicity. In the case of extrinsic s-enhanced stars, in 
addition to the initial masses and the metallicity of both the binary components, we should 
know the evolutionary status of the presently observed secondary star, its age at the epoch 
of the pollution and the total mass of chemically enriched material deposited on its surface 
{Bisterzo et al, submitted). Nevertheless, according to our low- mass- AGB models, the two 
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spectroscopic indexes [hs/ls] and [Pb/hs] should roughly attain the final value after a few 
third dredge-up episodes. In other words, if the s-process enhancement is sufficiently large, 
the knowledge of the precise evolutionary status is no longer a necessary condition. In 
practice, if the s-enhancement has been caused by the radiative ^'^C burning, as we have 
shown to occur in a low mass AGB stars, the [hs/ls] and the [Pb/hs] should only depend 
on the stellar metallicity and, eventually, on its mass. 

In the following analysis we consider s-enhanced stars spectroscopically classified as: 
S stars (0-rich AGB stars), N type C stars (C-rich AGB stars) post-AGB (either O- and 
C-rich), Ba and/or CH stars (C-rich extrinsic s-enhanced stars). We also consider a few SC 
stars, even though their evolutionary status and origin are unknown. Then, in Fig. [T2] we 
compare, as a function of the metallicity, our final surface [hs/ls] with the corresponding 
measured values, for a rather large sample of stars. In the caption we report the references 
corresponding to observational data. Although a similar trend with the metallicity may 
be recognized, the rather large spread of the observed [hs/ls] at any metallicity cannot be 
explained through a change of the stellar mass. 

The comprehension of the origin of this spread is a hard task. Note that the typical 
error associated to the derivation of the heavy elements abundances for these s-enhanced 
stars is of the order of ±0.1 dex, even if a more conservative evaluation might suggest up 
to ±0.3 dex. Such an error, however, only takes into account the uncertainty introduced 
in the abundance analysis by the unknown stellar parameters (effective temperature, 
gravity, metallicity). The systematic uncertainty associated to the particular physical 
prescriptions used to model the stellar atmosphere, such as the deviations from the local 
thermodynamics equilibrium (LTE) or from the spherical symmetry (ID versus 3D models), 
is often neglected. As a matter of fact, the use of different procedures for the abundance 
analysis and/or model atmospheres might lead to different stellar parameters and, hence, to 
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Fig. 12. — Our theoretical [hs/ls] curves compared with observational data. We re- 



port: intrinsic 0-stars (open blue pent agon: iSmith fc Lambert 



sic SC-stars (little o pen dark squares 



Abia &: Wallerstein 



filled dark squares: 



triangles: 



Abia et al. 



2002 



van Winckel fc Reyniers 



Zamora et al. 



20001 : 



Reddv et al. 



trinsic 0-rich stars (large filled gre en circ 



stars (little filled naagenta circles: 



Tomkin fc Lambert 



1986 



Smith 



es: 



1985 



1986 



19901 ): intrin- 



19981 ): intrinsic C-stars (little 



2009h: post-AGB stars (fille d red 



2002 



Revniers et al. 



Smith fc Lambert 



Smith fc Suntzeff 



984 



1985 



Smith fc Lambert 



1987 



Smith fc Lambert 



1986 



984 



2004); ex- 



1990f ) , Ba- 



Kovacs 



1990 



1995 



Pereira et al. 
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different results. A couple of examples may illustrate this problem. We found at least two 
stars in our sample for which a revision of the abundance analysis is possible. The SC-type 
star RR Her was analyzed by 



Abia Sz Wallerstein 



(Il998l ) by using atmosphere models 



from Alexander et al. (unpub 



[hs/ls]=0.5. Recently, 



Z amor a 



i shed) 



. For this star, they obtained [Fe/H]=0.2 and an 



20091) reanalyzed this s tar by using the state-of-the-art of 



atmosphere models for C-stars (iGustafsson et al. 



20081 ) ■ together with updated molecular 



and atomic line lists. The new derived [Fe/H] is -0.09 and the [hs/ls] ratio is 0.06. As a 
result, the revised [hs/ls] nicely fits the theoretical prediction (see Fig. [121). Moreover, w'e 
were also able to re-analyze the Ba-star HR 107, previously studied by I Allen &: Barbuy 



( I2OO6I ). We take advantage of a very high-resolution (R~ 200000) HARPS visual spectrum 



available in the 



Allen fc Barbuy 



SO ar chive. Then, by adopting the same atmosphere parameters as in 



(l2006l ). we derive a [hs/ls] =0.2 to be compared with the [hs/ls] =0 obtained 



by the auth ors. A more detai. 



be found in 



Abia et al. 



(12002), 



e d description of the tools used for the revised 



de Laverny et al. 



fl2006h and 



Zamora et al. 



analys is may 



fl2009h . 



Lead measurement s are also available for a sub-sample of the observed stars 



(^AUen fc Barbuy 



20061 ). 



The comparison with the theoretical [Pb/hs] index is shown 
in Fig. [131 As for the [hs/ls] ratio, the general trend of the observed data appears in 
reasonable agreement with the theoretical expectations. Also in this case, however, the 
observed spread cannot be ascribed to a variation of the stellar mass within the range 
1.5 — S.OMq. We have also revised the lead abundance of the Ba stars HR 107. From the 



HARPS high- resolution spectrum we have o btained an upper l imit 



is definitely smaller than the 0.34 quoted by 



-< r 



Allen fc Barbuy , (20061 ). This confirms the 



Pb/hs]=0.05, which 



need of more accurate evaluation of the systematic uncertainties affecting the abundance 
analysis. 



Although it is possible that the observed spread may be (totally or in part) ascribed 
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to a large observational uncertainty, we mention possible alternative theoretical scenarios, 
which may result in a certain spread of the spectroscopic indexes [hs/ls] and [Pb/hs]. As 
discussed in Section [31 these indexes basically depend on 3 quantities, namely: neutrons, 
poisons and seeds. A change of the amount of neutrons would imply a ch ange of the 



( effec tive) ^^C in the pocket. In the post-process calculations perforraed by 



1998 



fsee also 



Busso et al. 



2001 



Pignatari et al. 



2004 : 



Husti et al. 



2008 



Gallino et al. 



Bisterzo et al. 



201 ll ). the amount of ^^C within the pocket (and accordingly the amount of ^^N) has 
been varied parametrically in order to match the observed chemical pattern. In this way, 
since the extension of the pocket and, in turn, the amount of Fe seeds are unchanged, the 
neutrons-to-seed ratio can be freely varied, in order to reproduce the observed spread of the 
spectroscopic indexes. Such a procedure is justified by the uncertain theoretical description 
of the hydrodynamical process which determines the proton profile in the transition zone 
between the fully-convective envelope and the radiative H-exhausted core during a TDU 
episode. As a matter of fact, in our model such a profile is the result of the assumed 
exponential decay of the average mixing velocity at the inner border of the convective 
envelope. A spread of such an average description may eventually imply a spread in the 
resulting profile. Note, however, that even a small variation (positive or negative) of the P 
parameter determining the shape of the convective velocities within the convective-radiative 
transition region strongly depletes the resulting 5-process nucleosynthesis, as demonstrated 
in Section H] (see also Paper I). This would lead to definitely lower final surface heavy 
elements abundances, thus in contrast with observations. Moreover, the C-star luminosity 
function would be altered (see Section EH]). Such an occurrence provides an independent 
check of the adopted parametrization. As we have already shown, the predicted and the 
observed LF of the Galactic C-stars are in good agreement. 

Alternative scenarios may eventually account for the observed spread of the 
spectroscopic indexes, leaving the LFCS unchanged. Mixing induced by rotation has been 




Fig. 13. — Theoretical [Pb /hs] curves compared w ith the same ratios as observed in a sample 
of Galactic barium stars (lAUen fc Barbuy II2006I ). The new upper limit derived for HR107 
(this work) is also reported. 
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invoked to explain the formation of the ^^C pocket (ILanger et aL Ill999l ). The common 
conclusion of all these attempts is that the resulting ^^C pocket is too small for a sizeable 
s-process nucleosynthesis. This likely occurs because, during TDU episodes, the timescale 
for the formation of a proton profile (and the subsequent ^'^C pocket) in the upper layers of 
the He-intershell is too short compared to the timescale characterizing rotational-induced 
instabilities. Nevertheless, these instabilities might act during the longer time elapsed 

and the 



beginning of the s-process nucleosynthesis (iHerwig et al. 



e (maximum c 


epth 


(Herwie et al. 


2003 



Siess et al. 



20041 ). In 



low-mass AGB stars of nearly solar composition, such a period may be as long as a few 
10^ yr. If enough total angular momentum is conserved at the beginning of the TP-AGB 
phase, the Eddington- Sweet instability, in particular, may have enough time to smear the 
^^C (and the ^^N) pocket, thus reducing the local neutrons-to-seeds ratio. A quantitative 
analysis of this hypothesis is beyond the purpose of the present work and will be addressed 
in a forthcoming paper. 

Finally, let us note that when the s-process nucleosynthesis takes place at the base 
of a quite extended convective zone, the seeds consumption is counterbalanced by the 
rapid mixing of fresh material from outside. On the contrary, due to the rapid destruction 
compared to the production, the amount of neutrons is maintained at a local equilibrium 
value, which depends on the temperature. As a result, a relatively low neutrons-to-seed ratio 
is maintained compared to the case of radiative burning, thus favoring the accumulation of 
light-s elements. In massive AGB stars, where the main neutron source is the ^^Ne(a, n)^^Mg 
reaction, which is active at the base of the convective shell, we expect an overproduction of 
Is. The case of the 3 Mq model with Z = 0.001, where a substantial activation of the ^^Ne 
convective burning is found, shows in fact a particularly low [hs/ls] compared to lower mass 
models with the same metallicity. 
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7. Conclusions 

We presented FRUITY, an on-line database entirely dedicated to the nucle- 
osynthesis occurring in AGB stars. Its interface is hosted on the web-pages of 
the Teramo Observatory (INAF-OACTe) and is accessible at the internet address 
http: / /www.oa-teramo.inaf.it/fruity[ In this paper we discussed and analyzed the physical 
and chemical evolution of the first set of models included in the database, currently 
focused on low mass AGB stars (1.5 < M/Mq < 3.0) with metallicities in the range 
1 X 10~^ < Z < 2 X 10~^. In the future, we intend to extend the grid of models by 
increasing the mass (1.0 < M/Mq < 7.0) and metallicity (—3.2 < [Fe/H]< —1.2) ranges. 
Our calculations include a full nuclear network, from H to Bi, with 700 isotopes linked 
by more than 1000 nuclear processes. Thus, we followed in detail the nucleosynthesis of 
both light and heavy (s-process) elements. On the web interface, models selection can be 
operated choosing the desired mass (M) and metallicity {Z) values, while properties of 
nuclides (elements, isotopes, s-process spectroscopic indexes and yields) can be retrieved 
for all masses and atomic numbers or only for a sub-sample (selecting an appropriate pair 
of A/Z values). Moreover, the surface distributions after each TDU or just the final one 
can be selected for each model. Elemental abundances are given in the usual spectroscopic 
notation, while mass fractions are returned for the isotopes. Finally, net yields are given in 
solar mass units. 

We analyzed the effects deriving from the activation of both the neutron sources in 
our low mass AGB models, the ^^C(Q;,n)^^0 and the ^^Ne(a,n)^^Mg reactions, highlighting 
their contribution at different metallicities and for different initial masses. We evaluated 
the sensitivity of our model predictions to the adopted input parameters, with particular 
emphasis on the treatment of convection (the a parameter of the mixing length and the 
/3 parameter determining the shape of the exponentially decreasing velocity profile at the 
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inner border of the convective envelope) and the mass loss law. We compared our model 
with the observational counterparts, from both a physical (luminosity function of carbon 
stars) and chemical (s-process spectroscopic indexes at different metallicities) point of view. 
Our theoretical luminosity function of carbon stars well fits the observational one, recently 
derived with the state-of-the-art set of distances (HIPPARCOS) and spectra (ISO). If recent 
SPITZER data are taken into account, a good agreement was also found for the Magellanic 
Clouds. In the future this comparison could be improved for the C-star population of 
nearby galaxies, as soon as a larger number of IR data is available. A comparison between 
our theoretical curves and spectroscopic data showed that, at fixed metallicity, our models 
cannot account for the s-process indexes spread identified by observations (in particular the 
[hs/ls] ratio). In the framework of our current modelling, we excluded that such a spread 
could be ascribed to the different initial masses of AGB models. We speculated on possible 
solutions from both an observational and theoretical point of view. In this latter case, we 
suggested that rotation may (at least partially) solve the current disagreement. 

Part of this work was supported by the Spanish grants AYA2008-04211-C02-02 and 
FPA2008-03908 from the MEG and the Itahan Grant FIRB 2008 "FUTURO IN RIGERGA" 
G81J10000020001. We deeply thank R. Guandalini for providing unpubhshed data. We 
thank the anonymous referee for a detailed reading of this paper and for helpful comments 
and suggestions. 
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Table 1. The C-enhancement factors for the H-rich opacity tables at various metallicities, 
defined as the ratio between the ^^C abundance used to calculate the opacity table and the 

initial ^^C. 



Z [Fe/H] enhancement factor (/j) 
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Table 2. Fundamental properties of the pre-AGB phase. 
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Table 2 — Continued 
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Table 2 — Continued 
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Table 3. Evolution of selected physical parameters of the models at Z = 2 x 10 
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3. 


,18E- 


■02 


3 


.lOE- 


■03 


9, 


.31E- 


■01 


O.OOE+00 


12, 


.5 


2, 


.12 


2, 


.02E- 


■02 


0. 


,31 




-2W 


1, 


.89 


0, 


.549 


O.OOE+00 


3. 


,18E- 


•02 


5 


.40E- 


■03 


8, 


.66E- 


■01 


O.OOE+00 


16, 


.1 


2, 


.26 


2. 


,02E- 


■02 


0. 


,31 




-1 


1, 


.88 


0, 


.555 


O.OOE+00 


3. 


,12E- 


-02 


6 


.20E- 


■03 


8, 


.30E- 


■01 


O.OOE+00 


18, 


.0 


2, 


.39 


2, 


.02E- 


■02 


0. 


,31 




1 


1 


.86 


0, 


.562 


8.00E-04 


3. 


,00E- 


■02 


6 


.80E- 


■03 


7, 


.85E- 


■01 


1.18E-01 


18, 


.1 


2, 


.48 


2, 


.03E- 


■02 


0. 


,32 


17 


2 


1, 


.85 


0, 


.569 


1.70E-03 


2. 


,87E- 


-02 


8, 


.30E- 


■03 


7, 


.21E- 


■01 


2.05E-01 


17, 


.5 


2, 


.56 


2, 


.06E- 


■02 


0. 


,37 


33 


3 


1, 


.83 


0, 


.577 


3.10E-03 


2. 


,75E- 


-02 


9 


.lOE- 


■03 


6, 


.82E- 


■01 


3.41E-01 


16, 


.8 


2, 


.63 


2, 


.12E- 


■02 


0. 


,45 


36 


4 


1 


.80 


0, 


.583 


4.00E-03 


2. 


,65E- 


■02 


1 


.OOE- 


■02 


6, 


.32E- 


■01 


4.00E-01 


16, 


.1 


2, 


.70 


2, 


.20E- 


■02 


0. 


,57 


31 


5 


1, 


.76 


0, 


.590 


4.60E-03 


2. 


,55E- 


■02 


1 


.05E- 


■02 


5, 


.94E- 


■01 


4.38E-01 


15, 


.4 


2, 


.75 


2, 


.30E- 


■02 


0. 


,70 


29 


6 


1, 


.70 


0, 


.596 


4.90E-03 


2, 


,44E- 


-02 


1 


.08E- 


■02 


5, 


.69E- 


■01 


4.54E-01 


14, 


.5 


2, 


.79 


2, 


.41E- 


■02 


0, 


,84 


29 


7 


1, 


.61 


0, 


.602 


4.90E-03 


2. 


,33E- 


-02 


1 


.08E- 


■02 


5, 


.44E- 


■01 


4.54E-01 


13, 


.5 


2, 


.83 


2, 


.51E- 


■02 


0. 


,99 


28 


8 


1, 


.49 


0, 


.608 


4.70E-03 


2. 


,22E- 


■02 


1 


.07E- 


■02 


5, 


.25E- 


■01 


4.39E-01 


12, 


.4 


2, 


.85 


2, 


.63E- 


■02 


1. 


,15 


26 


9 


1, 


.34 


0, 


.614 


4.50E-03 


2. 


llE- 


-02 


1 


.04E- 


■02 


5, 


.17E- 


■01 


4.33E-01 


11, 


.4 


2, 


.87 


2, 


.76E- 


■02 


1. 


,33 


20 


10 


1, 


.17 


0, 


.619 


1.15E-03 


2. 


OOE- 


-02 


9 


.60E- 


■03 


5, 


.14E- 


■01 


1.20E-01 


10, 


.2 


2, 


.88 


2, 


.76E- 


■02 


1. 


,33 


19 



Table 3 — Continued 



n. 



TP 



AM^DU 



AM^ 



A 



At^ 



rpd 

-'-MAX 



-'surf 



C/0 



surf 



/l^Tot 
-'"TDU 



3.44E-02 Mf. 







Tc = 3.39E+05 yr 



M=2.5 M 





2.47 


0.523 


O.OOE+00 


3.76E-02 


6.70E-03 


8.52E-01 


O.OOE+00 


23.0 


2.25 


2.02E-02 


0.31 




1 


2.46 


0.531 


8.00E-04 


3.70E-02 


7.30E-03 


7.50E-01 


8.16E-02 


36.1 


2.42 


2.04E-02 


0.34 


5 


2 


2.45 


0.540 


1.50E-03 


3.39E-02 


9.00E-03 


7.48E-01 


2.05E-01 


25.8 


2.44 


2.05E-02 


0.36 


16 


3 


2.44 


0.548 


2.80E-03 


3.27E-02 


9.30E-03 


7.29E-01 


3.01E-01 


23.4 


2.52 


2.08E-02 


0.40 


30 


4 


2.43 


0.555 


4.00E-03 


3.15E-02 


l.OlE-02 


6.91E-01 


3.96E-01 


21.9 


2.59 


2.13E-02 


0.48 


35 


5 


2.41 


0.562 


5.20E-03 


3.04E-02 


1.09E-02 


6.53E-01 


4.77E-01 


20.7 


2.65 


2.20E-02 


0.57 


30 


6 


2.39 


0.568 


6.00E-03 


2.94E-02 


1.16E-02 


6.15E-01 


5.17E-01 


19.7 


2.70 


2.28E-02 


0.68 


27 


7 


2.37 


0.574 


6.60E-03 


2.84E-02 


1.21E-02 


5.86E-01 


5.45E-01 


18.6 


2.75 


2.37E-02 


0.80 


26 


8 


2.34 


0.580 


6.70E-03 


2.73E-02 


1.23E-02 


5.61E-01 


5.45E-01 


17.5 


2.78 


2.45E-02 


0.91 


24 


9 


2.30 


0.585 


6.50E-03 


2.61E-02 


1.22E-02 


5.43E-01 


5.33E-01 


16.3 


2.81 


2.53E-02 


1.02 


23 
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Ml, 




AM^Du 




AM^ 


r 


A 


Ai? 


rpd 

-'-MAX 




C/ Osurf 


ric 


10 


2.24 


0.591 


6.30E-03 


2.50E-02 


1.20E-02 


5.30E-01 


5.25E-01 


15.1 


2.84 


2.61E-02 


1.13 


22 


11 


2.16 


0.596 


6.10E-03 


2.38E-02 


1.18E-02 


5.16E-01 


5.17E-01 


14.0 


2.85 


2.69E-02 


1.24 


22 


12 


2.06 


0.602 


5.90E-03 


2.28E-02 


1.15E-02 


5.07E-01 


5.13E-01 


13.0 


2.87 


2.77E-02 


1.36 


22 


13 


1.92 


0.607 


5.70E-03 


2.19E-02 


1.12E-02 


4.99E-01 


4.09E-01 


12.0 


2.89 


2.85E-02 


1.48 


21 


14 


1.76 


0.612 


5.30E-03 


2.09E-02 


1.09E-02 


4.90E-01 


4.86E-01 


11.2 


2.91 


2.95E-02 


1.60 


21 


15 


1.57 


0.617 


4.90E-03 


2.00E-02 


1.04E-02 


4.89E-01 


4.71E-01 


10.3 


2.92 


3.05E-02 


1.75 


21 


16 


1.35 


0.623 


3.90E-03 


1.91E-02 


l.OOE-02 


4.86E-01 


3.90E-01 


9.5 


2.92 


3.15E-02 


1.89 


21 



M^^ty = 7.82E-02 Mq tc = 9.27E+05 yr 



M=3.0 M 



-3 2.98 0.550 O.OOE+00 2.96E-02 2.50E-03 9.54E-01 O.OOE+00 9.7 2.06 2.02E-02 0.31 
-2(*) 2.97 0.555 O.OOE+00 3.04E-02 4.40E-03 8.81E-01 O.OOE+00 13.4 2.24 2.02E-02 0.31 
-1 2.97 0.560 O.OOE+00 2.96E-02 5.50E-03 8.33E-01 O.OOE+00 15.3 2.35 2.02E-02 0.31 
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1 


2, 


.96 


0, 


.567 


8, 


.OOE- 


■04 


2, 


.89E- 


■02 


6. 


,80E-03 


7. 


,81E-01 


1 


.18E- 


■01 


16.0 


2. 


,46 


2, 


.02E-02 


0, 


.32 


2 


2 


2, 


.96 


0, 


.574 


2, 


.30E- 


■03 


2, 


.77E- 


■02 


7, 


,80E-03 


7. 


,34E-01 


2, 


.95E- 


■01 


15.9 


2. 


,55 


2, 


.04E-02 


0, 


.35 


16 


3 


2, 


.94 


0, 


.581 


3, 


.60E- 


■03 


2, 


.69E- 


■02 


9. 


,10E-03 


6. 


,74E-01 


3, 


.96E- 


■01 


15.7 


2. 


,63 


2, 


.08E-02 


0, 


.41 


35 


4 


2, 


.93 


0, 


.587 


4, 


.90E- 


■03 


2, 


.60E- 


■02 


9. 


,90E-03 


6. 


,31E-01 


4, 


.95E- 


■01 


15.4 


2. 


,71 


2, 


.14E-02 


0, 


.48 


33 


5 


2, 


.91 


0, 


.593 


5, 


.90E- 


■03 


2, 


.53E- 


■02 


1. 


,08E-02 


5. 


,81E-01 


5, 


.46E- 


■01 


15.0 


2. 


,77 


2, 


.20E-02 


0, 


.57 


31 


6 


2, 


.87 


0, 


.598 


6, 


.70E- 


■03 


2, 


.45E- 


■02 


1. 


,14E-02 


5. 


,45E-01 


5, 


.88E- 


■01 


14.5 


2, 


,82 


2, 


.26E-02 


0, 


.66 


27 


7 


2, 


.81 


0, 


.603 


7, 


.OOE- 


■03 


2, 


.37E- 


■02 


1. 


,18E-02 


5. 


,11E-01 


5, 


.93E- 


■01 


13.9 


2. 


,87 


2, 


.33E-02 


0, 


.76 


26 


8 


2, 


.74 


0, 


.608 


7, 


.OOE- 


■03 


2, 


.30E- 


■02 


1. 


,18E-02 


4. 


,91E-01 


5, 


.93E- 


■01 


13.3 


2. 


,91 


2, 


.40E-02 


0, 


.86 


25 


9 


2, 


.63 


0, 


.613 


7, 


.lOE- 


■03 


2, 


.21E- 


■02 


1. 


,17E-02 


4. 


,78E-01 


6, 


.07E- 


■01 


12.5 


2, 


,93 


2, 


.47E-02 


0, 


.96 


25 


10 


2, 


.49 


0, 


.617 


6, 


.80E- 


■03 


2, 


.13E- 


■02 


1. 


,16E-02 


4. 


.61E-01 


5, 


.86E- 


■01 


11.7 


2. 


,95 


2, 


.54E-02 


1, 


.07 


25 


11 


2, 


.32 


0, 


.622 


6, 


.50E- 


■03 


2, 


.05E- 


■02 


1. 


,13E-02 


4. 


,55E-01 


5, 


.75E- 


■01 


11.0 


2. 


,97 


2, 


.62E-02 


1 


.18 


24 


12 


2, 


.10 


0, 


.626 


6, 


.20E- 


■03 


1, 


.97E- 


■02 


1. 


,09E-02 


4. 


,52E-01 


5, 


.74E- 


■01 


10.3 


2. 


,99 


2, 


.71E-02 


1, 


.30 


22 


13 


1, 


.84 


0, 


.631 


5, 


.70E- 


■03 


1, 


.89E- 


■02 


1. 


,07E-02 


4. 


,41E-01 


5, 


.33E- 


■01 


9.6 


2, 


,99 


2, 


.80E-02 


1, 


.43 


22 


14 


1, 


.55 


0, 


.635 


4, 


.90E- 


■03 


1, 


.80E- 


■02 


1. 


,01E-02 


4. 


,44E-01 


4, 


.85E- 


■01 


8.8 


3. 


,00 


2, 


.91E-02 


1, 


.59 


21 



M^°* T = 7.54E-02 Mr 



Tc = 4.62E+05 yr 



Table 3 — Continued 



"Pulse number (negative values correspond to TP not followed by TDU) 
^10^ yr 

(*)Mh of this TP corresponds to M^^ of Tabled 



Table 4. Evolution of selected physical parameters of the models at Z = 6 x 10 



M=1.5 Mr 



-6 


1.6b 


rv ceo 

0.552 


O.OOb+OO 


3.13E-02 


O CATTl AO 

2.50E-03 


A C AT71 A 1 

9.59E-01 


A AATTl 1 AA 

O.OOE+00 


1 A O 

10.2 


O AT 

2.07 


G CiGTT^ AO 

6.06E-03 


A OO 

0.32 




-2W 


1.35 


0.556 


O.OOE+00 


3.18E-02 


4.50E-03 


8.81E-01 


O.OOE+00 


15.6 


2.26 


6.06E-03 


0.32 




-1 


1.34 


0.563 


O.OOE+00 


3.11E-02 


6.80E-03 


8.17E-01 


O.OOE+00 


17.9 


2.40 


6.06E-03 


0.32 




1 


1.33 


0.570 


2.00E-04 


2.97E-02 


6.90E-03 


7.87E-01 


2.90E-02 


18.1 


2.51 


6.06E-03 


0.32 


42 


2 


1.31 


0.578 


l.lOE-03 


2.78E-02 


8.30E-03 


7.29E-01 


1.33E-01 


17.2 


2.57 


6.25E-03 


0.42 


120 


3 


1.29 


0.585 


1.60E-03 


2.64E-02 


8.70E-03 


6.81E-01 


1.84E-01 


16.1 


2.65 


6.78E-03 


0.66 


130 


4 


1.26 


0.593 


2.20E-03 


2.50E-02 


9.20E-03 


6.39E-01 


2.39E-01 


14.9 


2.71 


7.62E-03 


1.05 


120 


5 


1.24 


0.600 


2.70E-03 


2.37E-02 


9.60E-03 


6.04E-01 


2.81E-01 


13.7 


2.75 


8.67E-03 


1.53 


110 


6 


1.20 


0.607 


2.80E-03 


2.24E-02 


9.60E-03 


5.78E-01 


2.92E-01 


12.5 


2.79 


9.84E-03 


2.05 


97 


7 


1.16 


0.614 


2.80E-03 


2.12E-02 


9.50E-03 


5.59E-01 


2.95E-01 


11.4 


2.82 


l.llE-02 


2.59 


91 


8 


1.09 


0.620 


2. lOE-03 


2.01E-02 


9.20E-03 


5.11E-01 


2.62E-01 


10.1 


2.81 


1.23E-02 


3.16 


87 



M;?'°* = 1.58E-02 M 



Tc = 7.86E+05 yr 
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M=2.0 M 



-2W 


1 


.97 


0, 


.555 


O.OOE+00 


3.32E- 


■02 


4.70E- 


■03 


8.82E-01 


O.OOE+00 


1 5.9 


2.27 


6.07E- 


■03 


0.29 




-1 


1 


.97 


0, 


.561 


O.OOE+00 


3.20E- 


•02 


6.30E- 


■03 


8.21E-01 


O.OOE+00 


2 0.3 


2.42 


6.07E- 


■03 


0.29 




1 


1 


.96 


0, 


.569 


1.60E-03 


3.06E- 


■02 


7.90E- 


■03 


7.58E-01 


2.03E-01 


2 0.9 


2.55 


6.29E- 


-03 


0.40 


110 


2 


1 


.95 


0, 


.577 


3.40E-03 


2.91E- 


■02 


9.50E- 


■03 


6.87E-01 


3.58E-01 


2 0.4 


2.65 


6.91E- 


■03 


0.69 


140 


3 


1 


.93 


0, 


.585 


5.30E-03 


2.80E- 


■02 


1.09E- 


■02 


6.25E-01 


4.86E-01 


1 9.6 


2.74 


7.91E- 


-03 


1.15 


110 


4 


1 


.91 


0, 


.591 


6.70E-03 


2.70E- 


■02 


1.20E- 


■02 


5.68E-01 


5.58E-01 


1 8.8 


2.82 


9.13E- 


-03 


1.70 


110 


5 


1 


.89 


0, 


.597 


7.30E-03 


2.60E- 


■02 


1.28E- 


■02 


5.19E-01 


5.70E-01 


17.7 


2.87 


1.04E- 


■02 


2.27 


94 


6 


1 


.86 


0, 


.603 


7.60E-03 


2.50E- 


■02 


1.29E- 


■02 


4.93E-01 


5.89E-01 


16.5 


2.90 


1.18E- 


■02 


2.85 


91 


7 


1 


.81 


0, 


.608 


7.60E-03 


2.40E- 


■02 


1.29E- 


■02 


4.71E-01 


5.89E-01 


15.3 


2.93 


1.31E- 


-02 


3.42 


87 


8 


1 


.75 


0, 


.613 


7.40E-03 


2.30E- 


■02 


1.27E- 


■02 


4.59E-01 


5.83E-01 


14.1 


2.96 


1.45E- 


-02 


3.99 


85 


9 


1 


.66 


0, 


.618 


7.10E-03 


2.20E- 


■02 


1.22E- 


■02 


4.58E-01 


5.82E-01 


12.9 


2.99 


1.58E- 


■02 


4.56 


84 
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Ml, 




AM^Du 




AM^ 


r 


A 


Ai? 


rpd 

-'-MAX 




C/ Osurf 


ric 


10 


1.55 


0.623 


6.50E-03 


2.10E-02 


1.18E-02 


4.45E-01 


5.51E-01 


11.9 


3.00 


1.72E-02 


5.13 


82 


11 


1.40 


0.628 


5.50E-03 


1.98E-02 


1.14E-02 


4.34E-01 


4.82E-01 


10.8 


3.00 


1.87E-02 


5.71 


80 


12 


1.21 


0.633 


4.50E-03 


1.87E-02 


1.06E-02 


4.45E-01 


4.25E-01 


9.68 


3.00 


2.02E-02 


6.35 


79 


13 


1.00 


0.638 


1.50E-03 


1.75E-02 


9.60E-03 


4.60E-01 


1.56E-01 


8.47 


2.99 


2.11E-02 


6.67 


65 


14 


0.83 


0.644 


7.00E-04 


1.54E-02 


7.80E-03 


5.04E-01 


9.00E-02 


6.91 


2.92 


2.17E-02 


6.91 


63 



^TDU = 7.27E-02 Mq tc = 1.38E+06 yr 



M=2.5 Mr 





2.47 


0.549 


O.OOE+00 


3.44E-02 


4.90E-03 


8.78E-01 


O.OOE+00 


18.1 


2.28 


6.06E-03 


0.29 




1 


2.47 


0.555 


3.00E-04 


3.30E-02 


6.70E-03 


8.16E-01 


4.48E-02 


22.1 


2.42 


6.06E-03 


0.29 


31 


2 


2.46 


0.563 


2.20E-03 


3.16E-02 


8.30E-03 


7.52E-01 


2.65E-01 


22.4 


2.54 


6.32E-03 


0.42 


140 


3 


2.45 


0.571 


4.50E-03 


3.04E-02 


9.90E-03 


6.85E-01 


4.55E-01 


21.9 


2.65 


6.94E-03 


0.72 


140 


4 


2.44 


0.578 


6.30E-03 


2.93E-02 


1.16E-02 


6.17E-01 


5.43E-01 


21.3 


2.74 


7.81E-03 


1.13 


110 
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5 


2, 


.42 


0, 


.585 


7, 


.60E- 


■03 


2, 


.84E- 


■02 


1 


.27E- 


■02 


5, 


.63E- 


■01 


5, 


.98E- 


■01 


20.5 


2, 


.81 


8, 


.82E-03 


1 


.60 


95 


6 


2, 


.40 


0, 


.591 


8, 


.50E- 


■03 


2, 


.75E- 


■02 


1, 


.35E- 


-02 


5, 


.20E- 


■01 


6, 


.30E- 


■01 


19.4 


2, 


.87 


9, 


.88E-03 


2, 


.08 


87 


7 


2, 


.38 


0, 


.596 


8, 


.80E- 


■03 


2, 


.66E- 


-02 


1, 


.39E- 


-02 


4, 


.86E- 


■01 


6, 


.33E- 


■01 


18.3 


2, 


.91 


1, 


.lOE-02 


2, 


.56 


75 


8 


2, 


.36 


0, 


.601 


9, 


.OOE- 


■03 


2, 


.56E- 


■02 


1 


.39E- 


■02 


4, 


.67E- 


■01 


6, 


.47E- 


■01 


17.1 


2, 


.95 


1 


.20E-02 


3 


.03 


73 


9 


2, 


.33 


0, 


.606 


9, 


.lOE- 


■03 


2, 


.47E- 


■02 


1, 


.38E- 


■02 


4, 


.49E- 


■01 


6, 


.59E- 


■01 


15.9 


2, 


.98 


1, 


.31E-02 


3 


.50 


73 


10 


2, 


.29 


0, 


.610 


8, 


.90E- 


■03 


2, 


.38E- 


-02 


1, 


.37E- 


-02 


4, 


.32E- 


■01 


6, 


.50E- 


■01 


14.8 


2, 


.99 


1, 


.41E-02 


3, 


.97 


72 


11 


2, 


.23 


0, 


.615 


8, 


.70E- 


■03 


2, 


.29E- 


-02 


1, 


.34E- 


-02 


4, 


.25E- 


■01 


6, 


.49E- 


■01 


13.8 


3 


.02 


1, 


.52E-02 


4, 


.43 


71 


12 


2, 


.16 


0, 


.619 


8, 


.50E- 


■03 


2, 


.20E- 


■02 


1 


.30E- 


■02 


4, 


.20E- 


■01 


6, 


.54E- 


■01 


12.8 


3 


.04 


1 


.63E-02 


4, 


.88 


69 


13 


2, 


.07 


0, 


.623 


8, 


.lOE- 


■03 


2, 


.12E- 


-02 


1, 


.27E- 


-02 


4, 


.llE- 


■01 


6, 


.38E- 


■01 


11.9 


3 


.05 


1, 


.73E-02 


5, 


.34 


66 


14 


1, 


.95 


0, 


.628 


7, 


.70E- 


■03 


2, 


.03E- 


-02 


1, 


.22E- 


-02 


4, 


.09E- 


■01 


6, 


.31E- 


■01 


11.1 


3, 


.06 


1, 


.85E-02 


5, 


.81 


79 


15 


1, 


.80 


0, 


.632 


7, 


.lOE- 


■03 


1 


.95E- 


■02 


1 


.18E- 


■02 


4, 


.05E- 


■01 


6, 


.02E- 


■01 


10.3 


3 


.06 


1 


.96E-02 


6, 


.29 


75 


16 


1, 


.62 


0, 


.636 


6, 


.30E- 


■03 


1, 


.87E- 


■02 


1, 


.12E- 


■02 


4, 


.lOE- 


■01 


5, 


.63E- 


■01 


9.5 


3 


.07 


2, 


.09E-02 


6, 


.80 


72 


17 


1, 


.39 


0, 


.640 


5, 


.lOE- 


■03 


1, 


.77E- 


-02 


1, 


.05E- 


-02 


4, 


.17E- 


■01 


4, 


.86E- 


■01 


8.7 


3 


.06 


2, 


.22E-02 


7, 


.35 


83 


18 


1, 


.13 


0, 


.644 


3 


.40E- 


■03 


1, 


.66E- 


-02 


9, 


.60E- 


-03 


4, 


.31E- 


■01 


3 


.54E- 


■01 


7.8 


3 


.04 


2, 


.36E-02 


7, 


.90 


70 


19 


0, 


.89 


0, 


.649 


1 


.lOE- 


■03 


1, 


.51E- 


■02 


8, 


.30E- 


■03 


4, 


.58E- 


■01 


1, 


.33E- 


■01 


6.5 


3 


.00 


2, 


.45E-02 


8, 


.25 


70 



Table 4 — Continued 



n. 



TP 



AM^DU 



AM^ 



A 



■'surf 



C/0 



surf 



/l^Tot 
-'"TDU 



1.21E-01 Mr: 







Tc = 1.99E+06 yr 



M=3.0 M 



-2 


2.96 


0.653 


O.OOE+00 


1, 


.67E- 


•02 










2.23 


6.06E- 


■03 


0.29 




.iW 


2.96 


0.656 


O.OOE+00 


1, 


.66E- 


■02 


3.40E-03 


8.24E-01 


O.OOE+00 


5.2 


2.41 


6.06E- 


■03 


0.29 




1 


2.95 


0.661 


8.00E-04 


1 


.63E- 


■02 


4.40E-03 


7.47E-01 


1.82E-01 


5.9 


2.57 


6.11E- 


■03 


0.31 


82 


2 


2.94 


0.665 


2.20E-03 


1 


.57E- 


■02 


5.40E-03 


6.74E-01 


4.07E-01 


6.1 


2.69 


6.31E- 


■03 


0.42 


100 


3 


2.93 


0.670 


3.20E-03 


1 


.54E- 


■02 


6.40E-03 


5.96E-01 


5.00E-01 


6.2 


2.79 


6.66E- 


■03 


0.59 


100 


4 


2.89 


0.673 


4.10E-03 


1 


.51E- 


■02 


7.00E-03 


5.44E-01 


5.86E-01 


6.3 


2.89 


7.13E- 


■03 


0.83 


100 


5 


2.84 


0.677 


4.90E-03 


1 


.49E- 


■02 


7.60E-03 


4.94E-01 


6.45E-01 


6.4 


2.97 


7.69E- 


■03 


1.10 


100 


6 


2.77 


0.680 


5.60E-03 


1 


.46E- 


■02 


8.20E-03 


4.50E-01 


6.83E-01 


6.4 


3.03 


8.32E- 


■03 


1.41 


99 


7 


2.66 


0.683 


5.90E-03 


1 


.44E- 


■02 


8.50E-03 


4.14E-01 


6.94E-01 


6.3 


3.09 


9.01E- 


■03 


1.75 


97 


8 


2.52 


0.686 


6.20E-03 


1 


.42E- 


■02 


8.70E-03 


3.92E-01 


7.13E-01 


6.2 


3.14 


9.77E- 


■03 


2.11 


94 
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9 


2.34 


0.689 


6.30E-03 


1.39E-02 


8.90E-03 


3.70E-01 


7.08E-01 


6.0 


3.16 


1.06E-02 


2.51 


82 


10 


2.11 


0.691 


6.10E-03 


1.36E-02 


8.80E-03 


3.59E-01 


6.93E-01 


5.8 


3.20 


1.16E-02 


2.97 


100 


11 


1.82 


0.694 


5.60E-03 


1.33E-02 


8.70E-03 


3.53E-01 


6.44E-01 


5.6 


3.21 


1.27E-02 


3.49 


110 


12 


1.45 


0.696 


4.10E-03 


1.28E-02 


8.30E-03 


3.58E-01 


4.94E-01 


5.3 


3.22 


1.39E-02 


4.05 


110 


13 


1.08 


0.700 


1.70E-03 


1.18E-02 


7.40E-03 


3.79E-01 


2.30E-01 


4.7 


3.20 


1.47E-02 


4.54 


100 



MT°tu = 5.67E-02 Mq tc = 4.93E+05 yr 



"Pulse number (negative values correspond to TP not followed by TDU) 
^10^ yr 

dlQ8 K 

(*)Mh of this TP corresponds to M^^ of Table [2] 
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Table 5. Evolution of selected physical parameters of the models at Z = 1 x 10 



M=1.5 Mr 



-2 


1 


.40 


0. 


,566 


O.OOE+00 


2. 


,81E- 


■02 


1. 


,70E- 


■03 


1.08E+00 


O.OOE+00 


7.8 


2. 


,00 


1 


.OlE-03 


0.29 




1 (*) 


1 


.69 


U, 


.biz 


A AATT^ 1 AA 

U.UUb+UU 


o 

6 


.loE- 


AO 

■02 


5 


.5UE- 


AO 


8.47E-01 


A AATT^ 1 AA 

O.OOb+UU 


17.0 


2. 


A A 

,40 


-I 

1, 


A 1 T7^ A O 

.Olb-03 


A OA 




1 


1 


.39 


0, 


.579 


2.30E-03 


3, 


.OlE- 


■02 


7, 


.50E- 


■03 


7.68E-01 


3.07E-01 


22.2 


2, 


.58 


1, 


.58E-03 


1.95 


740 


2 


1 


.38 


0, 


.587 


4.50E-03 


2, 


.84E- 


■02 


9, 


.70E- 


■03 


6.73E-01 


4.64E-01 


21.1 


2, 


.70 


3 


.OlE-03 


5.51 


830 


3 


1 


.36 


0, 


.593 


5.90E-03 


2, 


.71E- 


■02 


1, 


.12E- 


■02 


5.99E-01 


5.27E-01 


19.3 


2, 


.79 


4, 


.96E-03 


9.34 


710 


4 


1 


.34 


0, 


.599 


6.80E-03 


2, 


.59E- 


■02 


1, 


.21E- 


■02 


5.44E-01 


5.62E-01 


17.7 


2, 


.85 


7, 


.19E-03 


12.8 


680 


5 


1 


.31 


0, 


.605 


7.00E-03 


2, 


.48E- 


■02 


1 


.26E- 


■02 


5.03E-01 


5.63E-01 


16.2 


2, 


.90 


9, 


.33E-03 


15.5 


620 


6 


1 


.29 


0, 


.611 


7.10E-03 


2, 


.36E- 


■02 


1, 


.25E- 


■02 


4.86E-01 


5.60E-01 


14.6 


2, 


.93 


1 


.19E-02 


18.3 


610 


7 


1 


.25 


0, 


.616 


6.50E-03 


2, 


.24E- 


■02 


1, 


.21E- 


■02 


4.68E-01 


5.37E-01 


13.2 


2, 


.96 


1, 


.41E-02 


20.6 


590 


8 


1 


.22 


0, 


.621 


6.10E-03 


2, 


.13E- 


■02 


1, 


.16E- 


■02 


4.64E-01 


5.26E-01 


11.9 


2, 


.97 


1, 


.63E-02 


22.5 


560 


9 


1 


.19 


0, 


.626 


5.30E-03 


2, 


.OlE- 


■02 


1, 


.llE- 


■02 


4.56E-01 


4.77E-01 


10.7 


2, 


.98 


1 


.84E-02 


24.2 


510 


10 


1 


.15 


0, 


.631 


4.50E-03 


1 


.90E- 


■02 


1 


.04E- 


■02 


4.61E-01 


4.33E-01 


9.6 


2, 


.98 


2, 


.03E-02 


25.7 


530 



Table 5 — Continued 





Ml, 




AM^Du 




AM^ 


r 


A 


Ai? 


rpd 

-'-MAX 




C/ Osurf 


ric 


11 


1.11 


0.636 


3.60E-03 


1.78E-02 


9.80E-03 


4.62E-01 


3.67E-01 


8.6 


2.97 


2.20E-02 


26.9 


560 


12 


1.06 


0.642 


2.70E-03 


1.68E-02 


9.00E-03 


4.73E-01 


3.00E-01 


7.7 


2.96 


2.35E-02 


28.0 


520 


13 


1.00 


0.647 


1.80E-03 


1.57E-02 


8.30E-03 


4.79E-01 


2.17E-01 


6.9 


2.95 


2.47E-02 


28.8 


540 


14 


0.94 


0.653 


2.00E-04 


1.47E-02 


7.60E-03 


4.92E-01 


2.63E-02 


6.1 


2.93 


2.50E-02 


29.0 


480 



M^oty = 6.43E-02 Mq tc = 1.50E+06 yr 



M=2.0 Aff 



_l(*) 


1.94 


0.600 


O.OOE+00 


2.36E-02 


2.20E-03 


9.55E-01 


O.OOE+00 


6.48 


2.13 


l.OlE-03 


0.25 




1 


1.93 


0.604 


2.00E-04 


2.53E-02 


4.70E-03 


8.23E-01 


4.30E-02 


12.6 


2.47 


1.03E-03 


0.31 


370 


2 


1.92 


0.611 


3.00E-03 


2.41E-02 


6.70E-03 


7.31E-01 


4.48E-01 


15.1 


2.64 


1.56E-03 


1.85 


1000 


3 


1.91 


0.617 


5.40E-03 


2.33E-02 


8.80E-03 


6.27E-01 


6.14E-01 


15.0 


2.79 


2.63E-03 


4.56 


730 


4 


1.89 


0.622 


7.20E-03 


2.28E-02 


1.06E-02 


5.36E-01 


6.79E-01 


14.7 


2.90 


4.05E-03 


7.64 


670 


5 


1.88 


0.627 


8.50E-03 


2.23E-02 


1.19E-02 


4.72E-01 


7.14E-01 


14.2 


2.98 


5.64E-03 


10.6 


650 



Table 5 — Continued 





Ml, 


Mh 


AM^Du 






r 


A 


Af^ 


'-pd 

^ MAX 




C/ Osurf 




6 


1.85 


0.631 


9.30E-03 


2.19E-02 


1.26E-02 


4.27E-01 


7.38E-01 


13.6 


3.03 


7.31E-03 


13.3 


670 


7 


1.83 


0.635 


9.70E-03 


2.14E-02 


1.31E-02 


3.94E-01 


7.40E-01 


13.0 


3.08 


9.01E-03 


15.9 


650 


8 


1.81 


0.638 


9.70E-03 


2.09E-02 


1.31E-02 


3.78E-01 


7.40E-01 


12.2 


3.11 


1.07E-02 


18.3 


630 


9 


1.78 


0.641 


9.70E-03 


2.03E-02 


1.29E-02 


3.70E-01 


7.52E-01 


11.5 


3.15 


1.24E-02 


20.5 


630 


10 


1.75 


0.644 


9.50E-03 


1.97E-02 


1.27E-02 


3.58E-01 


7.48E-01 


10.8 


3.16 


1.41E-02 


22.5 


670 


11 


1.71 


0.647 


9.20E-03 


1.91E-02 


1.25E-02 


3.51E-01 


7.36E-01 


10.1 


3.18 


1.58E-02 


24.4 


640 


12 


1.66 


0.650 


8.80E-03 


1.85E-02 


1.20E-02 


3.51E-01 


7.33E-01 


9.46 


3.18 


1.74E-02 


26.2 


660 


13 


1.61 


0.653 


8.40E-03 


1.79E-02 


1.17E-02 


3.48E-01 


7.18E-01 


8.83 


3.18 


1.91E-02 


27.8 


650 


14 


1.53 


0.656 


8.00E-03 


1.72E-02 


1.13E-02 


3.50E-01 


7.08E-01 


8.24 


3.19 


2.08E-02 


29.4 


620 


15 


1.45 


0.659 


7.30E-03 


1.66E-02 


1.08E-02 


3.50E-01 


6.76E-01 


7.65 


3.18 


2.26E-02 


30.9 


640 


16 


1.34 


0.661 


6.40E-03 


1.58E-02 


1.02E-02 


3.56E-01 


6.27E-01 


7.06 


3.17 


2.44E-02 


32.3 


610 


17 


1.22 


0.665 


5.30E-03 


1.50E-02 


9.50E-03 


3.68E-01 


5.58E-01 


6.46 


3.17 


2.62E-02 


33.7 


600 


18 


1.09 


0.668 


3.70E-03 


1.41E-02 


8.70E-03 


3.82E-01 


4.25E-01 


5.83 


3.13 


2.80E-02 


34.9 


600 


19 


0.96 


0.672 


2.00E-03 


1.31E-02 


7.50E-03 


4.11E-01 


2.67E-01 


5.14 


3.10 


2.96E-02 


36.0 


610 



M^°\j = 1.31E-01 Mq tc = 1.85E+06 yr 



Table 5 — Continued 



M=2.5 M 



1 


2.45 


0.672 


5.00E-04 


1.56E-02 


3.10E-03 


8.30E-01 


1.61E-01 


4.8 


2.45 


1.02E-03 


0.26 


340 


2 


2.45 


0.676 


2.30E-03 


1.51E-02 


4.40E-03 


7.27E-01 


5.23E-01 


6.5 


2.65 


1.25E-03 


0.96 


440 


3 


2.44 


0.680 


4.20E-03 


1.49E-02 


6.00E-03 


6.14E-01 


7.00E-01 


7.0 


2.82 


1.81E-03 


2.52 


480 


4 


2.43 


0.683 


5.70E-03 


1.49E-02 


7.40E-03 


5.18E-01 


7.70E-01 


7.4 


2.96 


2.61E-03 


4.58 


600 


5 


2.42 


0.686 


6.80E-03 


1.49E-02 


8.50E-03 


4.44E-01 


8.00E-01 


7.7 


3.07 


3.56E-03 


6.78 


780 


6 


2.41 


0.688 


7.70E-03 


1.50E-02 


9.30E-03 


3.90E-01 


8.28E-01 


7.8 


3.16 


4.59E-03 


8.98 


820 


7 


2.38 


0.690 


8.20E-03 


1.50E-02 


9.90E-03 


3.44E-01 


8.28E-01 


7.8 


3.24 


5.67E-03 


11.1 


900 


8 


2.35 


0.693 


8.60E-03 


1.49E-02 


1.03E-02 


3.16E-01 


8.35E-01 


7.7 


3.30 


6.80E-03 


13.3 


910 


9 


2.30 


0.695 


9.00E-03 


1.49E-02 


1.06E-02 


2.96E-01 


8.49E-01 


7.6 


3.33 


7.97E-03 


15.3 


940 


10 


2.24 


0.696 


9.10E-03 


1.48E-02 


1.08E-02 


2.76E-01 


8.43E-01 


7.4 


3.37 


9.19E-03 


17.4 


1000 


11 


2.16 


0.698 


9.00E-03 


1.46E-02 


1.08E-02 


2.67E-01 


8.33E-01 


7.2 


3.40 


1.05E-02 


19.3 


1100 



Table 5 — Continued 



12 


2.07 


0.700 


9.10E-03 


1.45E-02 


1.07E-02 


2.67E-01 


8.50E-01 


7.0 


3.41 


1.18E-02 


21.3 


1100 


13 


1.95 


0.701 


8.90E-03 


1.42E-02 


1.07E-02 


2.55E-01 


8.32E-01 


6.7 


3.43 


1.32E-02 


23.2 


1100 


14 


1.81 


0.703 


8.60E-03 


1.40E-02 


1.05E-02 


2.55E-01 


8.19E-01 


6.5 


3.45 


1.47E-02 


25.1 


1100 


15 


1.64 


0.705 


8.30E-03 


1.37E-02 


1.02E-02 


2.59E-01 


8.14E-01 


6.2 


3.45 


1.64E-02 


27.2 


1100 


16 


1.43 


0.706 


7.50E-03 


1.33E-02 


9.90E-03 


2.56E-01 


7.58E-01 


5.8 


3.44 


1.83E-02 


29.3 


1100 


17 


1.18 


0.708 


5.80E-03 


1.26E-02 


9.30E-03 


2.70E-01 


6.24E-01 


5.3 


3.42 


2.07E-02 


31.8 


1100 


18 


0.94 


0.710 


2.40E-03 


1.14E-02 


7.90E-03 


3.10E-01 


3.04E-01 


4.5 


3.36 


2.29E-02 


33.8 


970 



M^^ty = 1.22E-01 Mo Tc = 1.05E+06 yr 



M=3.0 Mq 

-iW 2.95 0.777 O.OOE+00 7.35E-03 ••• ••• 2.35 l.OlE-03 0.27 

1 2.95 0.779 9.00E-04 7.49E-03 1.90E-03 7.85E-01 4.74E-01 1.8 2.55 1.06E-03 0.43 460 

2 2.94 0.780 1.90E-03 7.25E-03 2.50E-03 6.77E-01 7.60E-01 1.9 2.71 1.25E-03 0.97 740 



Table 5 — Continued 



3 


2, 


.94 


0, 


.782 


2, 


.70E- 


■03 


7. 


14E-03 


3. 


,10E-03 


5, 


.66E- 


■01 


8 


.71E- 


■01 


2. 


,1 


2, 


.85 


1 


.60E-03 


1.86 


720 


4 


2, 


.92 


0, 


.783 


3 


.30E- 


•03 


7. 


,18E-03 


3. 


,80E-03 


4, 


.87E- 


■01 


8 


.46E- 


■01 


2. 


,3 


2, 


.99 


2, 


.OOE-03 


2.86 


840 


5 


2, 


.89 


0, 


.784 


3, 


.60E- 


■03 


7. 


17E-03 


4. 


,30E-03 


4, 


.18E- 


■01 


8 


.37E- 


■01 


2. 


,4 


3, 


.10 


2, 


.43E-03 


3.90 


920 


6 


2, 


.85 


0, 


.785 


3 


.80E- 


■03 


7. 


18E-03 


4. 


,60E-03 


3 


.73E- 


■01 


8 


.26E- 


■01 


2. 


,5 


3 


.18 


2, 


.89E-03 


4.99 


1000 


7 


2, 


.79 


0, 


.786 


3 


.90E- 


■03 


7. 


24E-03 


4. 


,80E-03 


3 


.39E- 


■01 


8 


.13E- 


■01 


2. 


,5 


3 


.26 


3 


.39E-03 


6.13 


1200 


8 


2, 


.72 


0, 


.787 


4, 


.OOE- 


■03 


7. 


29E-03 


5. 


,00E-03 


3 


.19E- 


■01 


8 


.OOE- 


■01 


2, 


,6 


3, 


.32 


3, 


.94E-03 


7.32 


1300 


9 


2, 


.62 


0, 


.788 


4, 


.30E- 


■03 


7. 


35E-03 


5. 


,10E-03 


3 


.17E- 


■01 


8 


.43E- 


■01 


2. 


,6 


3, 


.38 


4, 


.53E-03 


8.57 


1600 


10 


2, 


.51 


0, 


.789 


4, 


.40E- 


■03 


7. 


,37E-03 


5. 


,30E-03 


2, 


.88E- 


■01 


8 


.30E- 


■01 


2. 


,6 


3 


.43 


5, 


.18E-03 


9.88 


1900 


11 


2, 


.38 


0, 


.790 


4, 


.50E- 


■03 


7. 


38E-03 


5, 


,40E-03 


2, 


.73E- 


■01 


8 


.33E- 


■01 


2. 


,6 


3, 


.46 


5, 


.90E-03 


11.3 


2300 


12 


2, 


.23 


0, 


.791 


4, 


.70E- 


■03 


7. 


,39E-03 


5. 


,50E-03 


2, 


.61E- 


■01 


8 


.55E- 


■01 


2. 


,6 


3, 


.49 


6, 


.73E-03 


12.7 


2900 


13 


2, 


.03 


0, 


.792 


4, 


.80E- 


■03 


7. 


,40E-03 


5. 


,70E-03 


2, 


.37E- 


■01 


8 


.42E- 


■01 


2. 


,6 


3 


.52 


7, 


.70E-03 


14.3 


3600 


14 


1, 


.79 


0, 


.793 


4, 


.70E- 


■03 


7. 


41E-03 


5. 


,70E-03 


2, 


.32E- 


■01 


8 


.25E- 


■01 


2. 


,5 


3 


.56 


8, 


.90E-03 


16.1 


5000 


15 


1, 


.51 


0, 


.794 


3, 


.90E- 


■03 


7. 


34E-03 


5, 


,60E-03 


2, 


.36E- 


■01 


6 


.96E- 


■01 


2, 


,5 


3 


.58 


1, 


.04E-02 


18.1 


5500 


16 


1, 


.21 


0, 


.795 


1, 


.90E- 


■03 


6. 


,95E-03 


5. 


,10E-03 


2, 


.78E- 


■01 


3 


.73E- 


■01 


2. 


,3 


3 


.54 


1, 


.13E-02 


19.7 


4700 



M^°\j = 5.73E-02 Mq tc = 1.44E+05 yr 



Table 5 — Continued 



"Pulse number (negative values correspond to TP not followed by TDU) 
^10^ yr 

(*)Mh of this TP corresponds to M^^ of Tabled 



Table 6. Final enrichments of the three s- process peaks ([Is/Fe], [hs/Fe] and [Pb/Fe]) and s- process indexes ([hs/ls] 
and [Pb/hs]) for different masses at selected metallicities. The final surface composition of some key light elements, 

such as the final surface C/0, C/N and ^^C/^^C ratios are given. 







Z = 2.0 X 10-2 






Z = 6.0 X 10-3 






Z = 


10-3 

J- KJ 




M/Mq 


1.5 


2.0 


2.5 


3.0 


1.5 


2.0 


2.5 


3.0 


1.5 


2.0 


2.5 


3.0 


[Is/Fe] 


0.40 


0.93 


1.12 


1.02 


0.89 


1.36 


1.44 


1.03 


1.26 


1.35 


1.21 


1.03 


[hs/Fe] 


0.10 


0.45 


0.58 


0.51 


1.40 


1.75 


1.79 


1.56 


1.91 


1.98 


1.72 


1.37 


[Pb/Fe] 


0.02 


0.19 


0.27 


0.24 


1.30 


1.59 


1.61 


1.46 


2.86 


2.87 


2.75 


2.62 


[hs/ls] 


-0.30 


-0.48 


-0.54 


-0.51 


0.50 


0.39 


0.34 


0.53 


0.65 


0.63 


0.52 


0.35 


[Pb/hs] 


-0.09 


-0.26 


-0.31 


-0.27 


-0.10 


-0.16 


-0.18 


-0.10 


0.95 


0.89 


1.03 


1.25 


[C/Fe] 


0.16 


0.42 


0.56 


0.48 


0.82 


1.17 


1.23 


0.94 


2.11 


2.18 


2.08 


1.78 


[N/Fe] 

L / 


0.32 


0.35 


0.39 


0.42 


0.36 


0.40 


0.44 


0.47 


0.63 


0.56 


0.50 


0.49 


[0/Fe] 


0.00 


0.00 


-0.02 


-0.03 


0.02 


0.03 


0.01 


-0.01 


0.35 


0.32 


0.25 


0.18 


[F/Fe] 


0.12 


0.40 


0.59 


0.50 


0.55 


1.06 


1.18 


0.81 


2.08 


2.22 


1.99 


1.51 


[Ne/Fe] 


0.05 


0.17 


0.27 


0.22 


0.20 


0.56 


0.67 


0.32 


1.42 


1.60 


1.38 


0.80 


[Na/Fe] 


0.05 


0.16 


0.28 


0.21 


0.11 


0.36 


0.46 


0.23 


1.13 


1.28 


1.01 


0.61 


[Mg/Fe] 


0.00 


0.00 


0.01 


0.01 


0.01 


0.07 


0.11 


0.05 


0.43 


0.72 


0.73 


0.42 


(C/0)surf 


0.71 


1.33 


1.89 


1.59 


3.15 


6.92 


8.25 


4.55 


29.0 


35.9 


33.9 


19.7 


(C/N)surf 


2.48 


4.22 


5.34 


4.08 


10.5 


21.3 


22.1 


10.9 


111 


151 


157 


70.7 



Table 6 — Continued 



Z = 2.0 X 10-2 Z = 6.0 X 10" 



Z = 10" 



(i2C/i3C)surf 50 98 141 115 239 575 693 349 3600 5193 4302 196 



00 



-88- 



Table 7. Yields (in M©) of selected isotopes for Z = 2.0 x 10"^ 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 M© M=3.0 M( 





-2.25E-02 


-6.08E-02 


-7.76E-02 


-9.01E-02 


^He 


2.04E-02 


4.79E-02 


6.06E-02 


7.36E-02 


12^ 


4.61E-04 


8.22E-03 


1.20E-02 


l.lOE-02 




5.38E-05 


1.12E-04 


9.81E-05 


1.27E-04 


14N 


9.86E-04 


2.66E-03 


2.91E-03 


4.16E-03 


15N 


-1.40E-06 


-4.19E-06 


-4.18E-06 


-5.53E-06 


160 


-1.88E-05 


-1.51E-04 


-6.59E-04 


-1.21E-03 




4.97E-06 


5.17E-05 


7.03E-05 


7.12E-05 


180 


-3.64E-06 


-1.14E-05 


-1.17E-05 


-1.54E-05 


19p 


1.03E-07 


1.45E-06 


2.40E-06 


2.25E-06 


20Ne 


-8.64E-07 


-4.24E-06 


-1.70E-06 


-2.85E-06 


2iNe 


2.58E-09 


8.55E-08 


3.61E-07 


7.49E-07 


22Ne 


1.21E-04 


1.34E-03 


2.21E-03 


2.03E-03 


23Na 


4.58E-06 


3.95E-05 


6.94E-05 


6.67E-05 




9.25E-07 


1.43E-05 


2.96E-05 


3.03E-05 


25 Mg 


-4.76E-07 


-3.63E-07 


5.18E-07 


3.56E-06 


26Mg 


4.66E-07 


5.91E-06 


1.05E-05 


1.40E-05 


26 Al 


1.76E-07 


7.92E-07 


7.87E-07 


7.27E-07 


3601 


2.96E-10 


3.07E-09 


3.28E-09 


4.78E-09 


^^Ca 


5.57E-10 


9.09E-09 


9.70E-09 


1.39E-08 


60pe 


6.74E-09 


5.26E-08 


8.23E-08 


4.33E-08 



-89- 



Table 7 — Continued 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 Mq M=3.0 M( 



89y 


1.29E-08 


1.90E-07 


2.84E-07 


2.81E-07 


107pd 


1.17E-10 


1.84E-09 


2.39E-09 


2.46E-09 




4.67E-10 


l.OlE-08 


1.46E-08 


1.45E-08 


205p]jj 


3.75E-12 


9.44E-11 


1.39E-10 


1.40E-10 


208 


1.42E-10 


7.23E-09 


1.07E-08 


1.21E-08 



-90- 



Table 8. Yields (in Mq) of selected isotopes for Z = 6.0 x 10 ^. 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 Mq M=3.0 M( 





-2.79E-02 


-6.59E-02 


-l.lOE-01 


-7.06E-02 


^He 


2.41E-02 


5.09E-02 


8.50E-02 


5.77E-02 


12^ 


2.85E-03 


1.26E-02 


2.09E-02 


1.07E-02 




1.49E-05 


2.01E-05 


2.64E-05 


3.46E-05 


14N 


3.34E-04 


5.98E-04 


1.03E-03 


1.41E-03 


15N 


-4.55E-07 


-9.14E-07 


-1.37E-06 


-1.71E-06 


160 


6.01E-05 


1.69E-04 


6.84E-05 


-2.35E-04 




3.61E-06 


2.31E-05 


2.58E-05 


2.62E-05 


180 


-1.32E-06 


-2.75E-06 


-4.21E-06 


-5.05E-06 


19p 


2.39E-07 


1.59E-06 


3.05E-06 


1.43E-06 


20Ne 


-1.16E-06 


1.28E-06 


8.31E-06 


-3.67E-06 


2iNe 


2.96E-08 


2.54E-07 


7.32E-07 


5.70E-07 


22Ne 


1.67E-04 


1.20E-03 


2.33E-03 


7.65E-04 


23Na 


3.06E-06 


1.75E-05 


3.60E-05 


2.02E-05 




2.99E-06 


2.37E-05 


5.32E-05 


9.17E-06 


25 Mg 


1.85E-07 


6.12E-06 


1.58E-05 


1.28E-05 


26Mg 


5.28E-07 


7.02E-06 


1.75E-05 


1.79E-05 


26 Al 


8.57E-08 


1.23E-07 


1.35E-07 


2.31E-07 


3601 


2.59E-10 


1.41E-09 


2.17E-09 


2.59E-09 


^^Ca 


6.65E-10 


3.36E-09 


4.96E-09 


6.24E-09 


60pe 


1.69E-08 


5.00E-08 


2.28E-07 


4.31E-08 



-91 - 



Table 8 — Continued 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 Mq M=3.0 M( 



89y 


1.67E-08 


l.OOE-07 


1.74E-07 


7.13E-08 


107pd 


3.16E-10 


1.68E-09 


2.88E-09 


1.26E-09 




1.17E-08 


5.38E-08 


8.62E-08 


5.11E-08 


205p]jj 


2.96E-10 


1.04E-09 


1.45E-09 


7.54E-10 


208 


3.37E-08 


1.50E-07 


2.33E-07 


1.69E-07 



-92 - 



Table 9. Yields (in Mq) of selected isotopes for Z = 10 ^. 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 M© M=3.0 M( 





-6.11E-02 


-1.15E-01 


-1.05E-01 


-6.03E-02 


^He 


4.82E-02 


8.85E-02 


8.00E-02 


4.57E-02 


12^ 


l.llE-02 


2.23E-02 


2.20E-02 


1.30E-02 




2.86E-06 


3.63E-06 


4.75E-06 


9.73E-05 


14N 


9.81E-05 


1.52E-04 


1.83E-04 


2.31E-04 


15N 


-7.51E-08 


-1.65E-07 


-2.44E-07 


-2.84E-07 


160 


2.91E-04 


4.47E-04 


3.89E-04 


3.05E-04 




1.71E-06 


7.09E-06 


1.07E-05 


5.29E-06 


180 


-2.67E-07 


-5.07E-07 


-7.67E-07 


-8.26E-07 


19p 


1.71E-06 


3.97E-06 


2.97E-06 


1.60E-06 


20Ne 


6.35E-06 


2.80E-05 


1.57E-05 


3.05E-06 


2iNe 


2.45E-07 


l.OOE-06 


1.22E-06 


4.36E-07 


22Ne 


1.05E-03 


2.69E-03 


1.87E-03 


5.91E-04 


23Na 


1.33E-05 


3.25E-05 


2.24E-05 


1.03E-05 




2.19E-05 


7.33E-05 


2.73E-05 


3.35E-06 


25 Mg 


6.95E-06 


4.36E-05 


6.29E-05 


3.21E-05 


26Mg 


6.62E-06 


3.14E-05 


9.37E-05 


4.81E-05 


26 Al 


2.57E-08 


5.40E-08 


6.74E-08 


2.86E-08 


3601 


1.91E-10 


4.34E-10 


5.05E-10 


4.54E-10 


^^Ca 


3.77E-10 


8.58E-10 


9.62E-10 


6.83E-10 


60pe 


2.33E-08 


1.27E-07 


5.40E-07 


5.93E-07 



-93- 



Table 9 — Continued 



Isotope M=1.5 Mq M=2.0 Mq M=2.5 Mq M=3.0 M( 



89y 


7.08E-09 


1.55E-08 


1.48E-08 


l.OOE-08 


107pd 


l.llE-10 


1.99E-10 


1.66E-10 


1.25E-10 




5.50E-09 


1.08E-08 


7.89E-09 


4.39E-09 


205p]jj 


2.36E-10 


3.95E-10 


1.98E-10 


1.42E-10 


208 


4.49E-07 


8.48E-07 


8.66E-07 


7.95E-07 



Table 10. Surface s-process branched number isotopic ratios for different initial masses at selected metallicities. The 

tabulated ratios relate to the branchings occurring at '^'^Kr, '^'^Rb, ^^Zr, ^^^Xe and ^"^^Ce, respectively. Solar ratios are 

also reported. 



Isot. Ratio 


© 


Z = 2.0 X 10- 


■2 


Z = 6.0 X 10" 


3 




Z = 


10-^ 




M/Mq 




1.5 


2.0 


2.5 


3.0 


1.5 


2.0 


2.5 


3.0 


1.5 


2.0 


2.5 


3.0 


86Ki-/84Kr 


0.30 


0.27 


0.21 


0.18 


0.20 


0.37 


0.36 


0.43 


0.54 


0.46 


0.66 


0.91 


1.09 


«^Rb/«^Rb 


0.38 


0.38 


0.30 


0.27 


0.30 


0.53 


0.57 


0.85 


0.89 


0.81 


1.29 


1.61 


1.79 


962r/94zr 


0.16 


0.08 


0.03 


0.02 


0.03 


0.10 


0.07 


0.13 


0.17 


0.19 


0.43 


0.84 


1.18 


i34xe/i32xe 


0.36 


0.29 


0.14 


0.11 


0.13 


0.08 


0.06 


0.09 


0.18 


0.15 


0.46 


0.90 


0.79 


142Cg/140Ce 


0.13 


0.10 


0.04 


0.03 


0.03 


0.01 


0.01 


0.02 


0.05 


0.04 


0.14 


0.28 


0.30 



Table 11. Selected physical and chemical quantities for a model with M=2 Mq and Z = 6 x 10~^ (FRANEC 
Standard Model, FSM). Normalized percentage variations as a function of different input parameters arc reported. 
Percentage values are defined as the difference between Test cases and Standard Case, normalized to the Standard 

one: (Test-FSM)/SM*100. 







M?, 




Yield" f ^2 


Yield" f^^Y) 


Yield" f^^^La) 


Yield" f2°^Pb) 


fhs/lsl 


FSM 


0.644 


-4.98 


1.38E+06 


1.26E-02 


l.OOE-07 


5.38E-08 


1.50E-07 


0.39 




+0.5% 


-13% 


+18% 


+39 % 


+54% 


+38% 


+31% 


-15% 




-1.6% 


+17% 


-22% 


-36% 


-47% 


-40% 


-42% 


+13% 


MAGb 


-3.1% 


+20% 


-28% 


-30% 


-29% 


-24% 


-21% 


7% 


MAGB 


+2.1% 


-17% 


+21% 


32% 


30% 


24% 


17% 


-5% 


Pi 


-1.2% 


+2% 


6% 


16% 


-66% 


-79% 


-93% 


-53% 


^2 


+1.0% 


-9% 


-16% 


-19% 


-84% 


-87% 


-93% 


-36% 



" In M(7) units 



-96- 



-97- 



Table 12. Initial abundances and net stellar yields relative to a 3 Mq model with Z=0.02 



taken from our database (FRUITY) and from 



Karakas 



20101 (KIO). Percentage differences 



are also shown. 



Isotope Yields" (FRUITY) Yields" (KIO) [(FRUITY/K10)-1]*100 





-9 OlE-02 


-6 95E-02 


-29 


^He 


7 36E-02 


5 13E-02 


+43 


12Q 


1 lOE-02 


1 20E-02 


-9 


13Q 


1.27E-04 


l.lOE-04 


+16 


14N 


4 16E-03 


3 09E-03 


+35 


15N 


-5 53E-06 


-5 03E-06 


-10 


160 


-1 21E-03 


-9 02E-04 


-34 




7 12E-05 


4 68E-05 


+52 


ISQ 


-1 54E-05 


-1 37E-05 


-12 


19p 


2 25E-06 


4 13E-06 


-46 


on TV T 


-2.85E-06 


l.OlE-05 


-128 


2iNe 


7.49E-07 


1.05E-06 


-29 


22Ne 


2.03E-03 


3.45E-03 


-41 


23Na 


6.67E-05 


9.28E-05 


-28 


24Mg 


3.03E-05 


1.36E-05 


+123 


25Mg 


3.56E-06 


2.17E-05 


-84 


26 Mg 


1.40E-05 


2.58E-05 


-46 


27A1 


2.58E-06 


2.58E-06 


-19 


28Si 


-2.78E-07 


-2.78E-07 


-212 


29Si 


2.96E-07 


2.96E-07 


-19 


30Si 


1.28E-06 


1.28E-06 


+10 



Table 12— Continued 



Isotope Yields'^ (FRUITY) Yields'^ (KIO) [(FRUITY/K10)-1]*100 



P 8.06E-07 8.06E-07 -50 



" In M(7) units 



